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ABSTRACT 



Context. We present maps of the main cooling lines of the neutral atomic gas ([OI] at 63 and 145 jim and [CII] at 158 jim) and in the 
[OIII] 88 iim line of the starburst galaxy M82, carried out with the PACS spectrometer on board the Herschel satellite. 
Aims. Our aim is to study the nature of the neutral atomic gas of M82 and to compare this gas with the molecular and ionized gas in 
the M82 disk and outflow. 

Methods. By applying PDR modeling we were able to derive maps of the main ISM physical parameters, including the optical depth 
(r [C //]), at unprecedented spatial resolution (~ 300 pc). 

Results. We can clearly kinematically separate the disk from the outflow in all lines. The t [C ii] is less than 1 everywhere, is lower in 
the outflow than in the disk, and within the disk is lower in the starburst region. 

The [CII] and [OI] distributions are consistent with PDR emission both in the disk and in the outflow. Surprisingly, in the outflow, 
the atomic and the ionized gas traced by the [OIII] line both have a deprojected velocity of ~ 75 km s _1 , very similar to the average 
velocity of the outflowing cold molecular gas (~ 100 km s _1 ) and several times smaller than the outflowing material detected in Ha 
(~ 600 km s _1 ). This suggests that the cold molecular and neutral atomic gas and the ionized gas traced by the [OIII] 88 iim line are 
dynamically coupled to each other but decoupled from the Ha emitting gas. 

Conclusions. We propose a scenario where cold clouds from the disk are entrained into the outflow by the winds where they likely 
evaporate, surviving as small, fairly dense cloudlets (n H ~ 500 - 1000 cm -3 , G ~ 500 - 1000, T gas ~ 300 K). We show that the 
UV photons provided by the starburst are sufficient to excite the PDR shells around the molecular cores and probably also the ionized 
gas that flows at the same PDR velocity. The mass of the neutral atomic component in the outflow is > 5-12xl0 7 M© to be 
compared with that of the molecular component (3.3 x 10 8 M©) and of the Ha emitting gas (5.8 x 10 6 M©). The mass loading factor, 
Moutflow/SFR, of the molecular plus neutral atomic gas in the outflow is ~ 2. Energy and momentum driven outflow models can 
explain the data equally well, if all the outflowing gas components are taken into account. 

Key words, galaxies: individual (M82)- Infrared: ISM - Galaxies: ISM - Galaxies: kinematics and dynamics- Galaxies: starburst - 
Galaxies: starburst 



1. Introduction 



In the last decades, the primary role of galactic outflows on 
galaxy evolution has become evident (Heckman, Armus and 
Mileydgg Ol Mc Keith et al. lT995l Lehnert and Heckman [19951 
Heckman 1998, Veilleux, Cecil and B land-Hawthorn 120051 
Springel, Di Matteo and Hernquist 12005b . Outflows in galaxies 
powered by either Active Galactic Nuclei (AGNs) or starbursts 
have been introduced in theoretical models (e.g. Croton et al. 
2006 ) as sources of negative feedback, that terminate black-hole 
growth and/or star-formation and to explain the tight correlation 
between the mass of the black hole, and the velocity dispersion 
or mass of the bulge in which it resides (Magorrian et al. 1998, 
Gebhardt et al. 120001 Merritt and Ferrarese l200Tl Tremaine et al. 
2002, Giiltekin et al. [20091) . These outflows contain a hot (10 7 
K) metal-enriched wind, in addition to entrained cooler gas and 
dust. Therefore, they may significantly influence the chemical 
evolution of both the interstellar medium (ISM) of the galaxy 



and the intergalactic medium (IGM) where this processed mate- 
rial might be eventually dispersed. 

Nearly all galaxies with high star formation r ates sh ow signs 
of outflows (Rupke, Veilleux and Sanders l2005al 12005b l[20Qgcb . 
They also seem to be ubiquitous in Lyman Break galaxies at z 
~ 3 where they are capable of creating holes of ~ 100 kpc in 
the surroundi ng IG M (Pettini et al. 2001 , Shapley et al. [20031 
Steidel et al. 120101 Genzel et al. 2010), in Lyman a emitting 
galaxies at z ~ 4 (Finkelstein et al. 2011) and compact gas rich 
Damped Lyman-a at z= 2.2 (Noterdaeme et al. 2012 ). Recently, 
very massive molecular and neutral gas AGN-driven outflows 
have been discovered in Mrk 231 (Fischer et al. 2010, Feruglio 
et al. l2PT01 . Rupke and Veilleux l20TTl Cicone et al.[20T2t , in the 
local early type ga laxy N GC 1266 (Aalto et al. l2PTTi in a z = 2 
QSO (WeiB et al. 120121) and in others Ultra Luminous Infrared 
Galaxies (U LIRGs ) (Sturm et aLlSmJand Chung et al. [20TTT) . 
Sturm et al. (120111) provided the first observational evidence that 
molecular outflows in some ULIRGs are powered by the AGN 
and that the mass outflow rates of the outflowing material exceed 
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the Star Formation Rate (SFR) of the host galaxies, confirming 
the theoretically postulated fundamental role of the outflows on 
galaxy evolution. 

The evident cosmological importance of this phenomenon 
has brought many astronomers to study a few nearby galactic 
winds in great detail. M82 is the best studied case. It is the clos- 
est (at a distance of 3.6 Mpc, Freedman et al. 119941) and the 
brightest galaxy hosting a spectacular bipolar outflow along its 
minor axis, although not as massive and energetic as those re- 
cently discovered in ULIRGs. Nevertheless, thanks to its prox- 
imity and to its favorable disk inclination (i = 81°), M82 is one 
of the nearby galaxies most studied at all wavelengths and has 
therefore become the archetype of a starburst with an outflow, 
with no AGN. 

M82 has undergone two powerful starburst episodes located 
in the central ~ 500 pc of the galaxy. The first starburst was trig- 
gered 10 7 yr ago by the encounter with M81, 10 8 yr ago. The sec- 
ond starburst was most likely bar driven 5xl0 6 yr ago (Forster 
Schreiber et al. 2003 ). Both of these starbursts comprise sev- 
eral star clusters each hosting many hundreds of young massive 
stars (Melo et al. 12005b . The superwind which feeds the bipo- 
lar outflow along the minor axis of the galaxy is generated in 
the inner 350 pc of the galaxy where intense diffuse hard X-ray 
emission is detected (Strickland and Heckman 2007 ). The out- 
flow reveals filamentary structures at many wavelengths: they 
are well visible in any high spatial resolution observation such 
as in the optical Ha, [Nil] and [OIII] lines (Shopbefland Bland- 
Hawthorn 1998), Xray (Strickland and Heckman 2007 ), warm 
molecular hydrogen (Veilleux, Rupke and S waters 2009) and 
PAHs (Engelbracht et al. 2006 ). Cold and warm dust and molec- 
ular gas are found in a biconical-like structure extending up to 3 
kpc away from the stellar disk on both sides, as revealed by UV 
scattered light (Hoopes et al. 2005 ), PAH emission (Engelbracht 
et al.12006, Kaneda et al. 2010 ) sub-millimeter continuum emis- 
sion (Leeuw and Robson [20091 . Roussel et al. l20T0l), 12 CQ (/ 
l->0) (Walter et al. l2002l Wei B, Walt er and Scoville l2055l) . 12 CO 
(/ 3— >2) (Seaquist and Clark 2001 ) line emission and H2 near 
infrared observations (Veilleux, Rupke and S waters 2009 ). Most 
likely, this material has been entrained by the wind. Further out, 
most of the dust (Roussel et al. [20T0l) and gas (Yun et al. [T994t 
associated with the halo seems to originate from the tidal inter- 
action with M8 1 . 

The bibolar outflow shows different morphologies in the north 
and in the south, especially at optical wavelengths (Shopbell 
and Bland-Hawthorn, H9981 Westmoquette et al. I2009al l2509bl 
Seaquist and Clark l200"Tb . The southern wind is approaching us 
and is more clearly visible because the inclination of M82's disk 
provides a more direct sight-line to the southern part of the disk. 
The northern outflow is receding from us and is likely to be more 
heavily obscured by the disk itself. This is why most detailed 
studies of M82's outflow have focused on the southern rather 
than th e northern outflow of M82 (Westmoquette et el. 2009a, 
I2009bl Shopbell and Bland-Hawthorn [19981). 
The most accepted interpretation of the outflowing material in 
M82 is the following: the supernova explosions and stellar winds 
of the starburst blow away a collimated hot (~ 10 8 K) gas, de- 
tected in the center of the galaxy in X-ray (Griffiths et al. 120001) . 
The wind also drags cold material from the disk into the out- 
flow. The collision between the hot fluid and the entrained cold 
clouds produces X-ray emission further out along the outflow. 
These clouds are heated by shocks and the UV radiation from 
the starburst producing most of the Ha emission. The cold ma- 
terial in the winds has been observed in its dust and molecular 
component and is directly illuminated by the light of the star- 



burst for at least ~ 1 kpc in both directions (Veilleux, Rupke and 
Swaters 120091) . 

So far little observations are available on the neutral atomic 
component of the outflow. Observations in the HI line of the 
M81/M82 group revealed many tidal streams from M82 (Yun 
et al. [T9941 Chynoweth et al. l2008i but no detailed study has 
been conducted on the HI associated with the winds. Colbert 
et al. (1999) have studied the far infrared (FIR) fine struc- 
ture lines observed with the LWS instrument on board of the 
Infrared Space Observatory (ISO) but the spatial (~ 80") 
and spectral (~ 10 3 km s" 1 ) resolution were too poor to distin- 
guish the outflow and the disk components. This is why we ob- 
served M82 with the PACS spectrometer (Poglitsc h et al . 120101) 
on board of the Herschel satellite (Pilbratt et al. 2010 ) in the 
[OI] lines at 63.2 and 145.6 /zm, and the [Nil], [CII] and [OIII] 
lines at 121.89, 157.7 and 88.3 yum respectively, as part of the 
SHINING Guaranteed Time Key Project (P.I. E. Sturm). With 
these observations we were able to resolve for the first time both 
spatially and kinematically the disk and the outflow of M82 in 
the main FIR fine structure lines. 

This paper is organized as follows: details on the observations 
and data reduction are given in Section 2. The resulting maps and 
related uncertainties are presented in Section 3. In Section 4 we 
explain the Photo-dissociated Regions (PDR) modeling of our 
observations and we discuss the meaning of the physical param- 
eters resulting from this modeling. We also compare the various 
gas phases participating to the outflow and we discuss their en- 
ergetics and their relation. In Section 5 we summarize the main 
conclusions. 

2. PACS spectrometer observations and data 
reduction 

We have observed M82 with a 5x5 (4x4) raster in the blue 
(red) channel, and with step sizes in YZ spacecraft coordinates 
equal to: step z = 16.0" (24.0") and step Y 14.5" (22.0 "), in 
bright line scan chopped mode, with large chopper throw in the 
following fine structure lines: [OI] ( 3 P 2 ^ 3 Pi) and ( 3 Pi-> 3 P ) 
lines at 63.2 and 145.6 /mi, [CII] ( 2 Pi/ 2 ^ 3 P3/2) line at 157.7 fim 
and [OIII] ( 3 Po— > 3 Pi) line at 88.3 jim. The on source integration 
times per raster position are: ~ 40, 40, 65, and 33 s for the 
[OI] 63, [OIII] 88, [OI] 145 and [CII] 158 lines respectively. 
Althought the bright line mode executes less grating steps 
per scan than the faint mode, it nevertheless attains Nyquist 
sampling. The spectral resolutions for the lines are: ~ 3300 (87 
km s" 1 ), ~ 2400 (125 km s" 1 ), 1160 (257km s" 1 ) and ~ 1255 
(240 km s" 1 ), respectively. 

The spatial steps used in the rasters, in combination with the 
spaxel pattern of the instrument, ensure spatially Nyquist sam- 
pled final maps. 

For each line, each raster position has been reduced in HIPEQ 
with the standard pipeline up to the rebinnedCube task. Then, a 
polynomial continuum plus a Gaussian profile in each pixel at 
each raster position have been fitted, and the corresponding in- 
tensity, wavelength peak, Full Width Half Maximum (FWHM) 
and celestial coordinates were recorded. In an external package 
we built the final maps applying the drizzle algorithm (Fruchter 
and Hook 2002 ). Thus, we produced intensity, peak velocity and 
line width maps for the four observed lines. Moreover, thanks to 

1 HIPE is a joint development by the Herschel Science Ground 
Segment Consortium, consisting of ESA, the NASA Herschel Science 
Center, and the HIFI, PACS and SPIRE consortia. 
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Fig. 1. Left column: for each line the signal to noise ratios maps. Right column: for each line the "error" maps in W m 2 , obtained 
as explained in Section 2 with the central region masked. 



the brightness of M82 we were also able to produce high spa- 
tial resolution continuum maps using the continuum measured 
from our line observations as well as from the data in the paral- 
lel channels. The final maps cover the disk and the base of the 
outflow of M82 up to 1 kpc away from the disk in both direction, 
at an unprecendent spatial resolution (from 130 to 270 pc) and 
sensitivity (see Figure 12 of Poglitsch et al. 2010 ) at these wave- 
lengths. Taking into account pointing offset and jittering during 
observations and the uncertainties related to the fact that we ar- 
bitrarily assign the coordinates at the center rather than at the 



corner of the final pixel, we estimate an overall astrometric un- 
certainty of ~ 4 ". 

In addition to the rasters, we have observed M82 also in the 
( 3 Pi^ 3 P 2 ) [Nil] line at 121.89 jim in faint chop-nod line 
scan mode, with 3 separate staring observations pointed on the 
galaxy's center, the northern and the southern outflow. These 
data were reduced in HIPE up to the rebinned cube step resulting 
in 25 line intensity values per pointing. 

For each of the fitted spectra we have estimated the errors 
from the residual spectrum, i.e. the original spectrum minus the 
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Fig. 2. Integrated line emission of the 4 observed lines at their original spatial resolution. Typical PSF widths are shown for 
the shortest and longest wavelengths. Units are in Wm~ 2 . The black cross corresponds to the adopted center of M82, i.e. 
9h55m52.2s 69d40m46.6s (J2000) (Dietz et al [T986l) . North is up and East is towards the left. 



fitted continuum plus Gaussian profile. These errors contain the 
noise of the spectrum as well as systematic deviations from a 
perfect fit. These errors were quoted in the same units as the 
fitted line fluxes, i.e., W m~ 2 per each spaxel. The line flux er- 
ror maps were generated with exactly the same drizzle scheme 
as for the line map. In order to create the error maps by the 
same method as for the flux, we used the intermediate quan- 
tity 1/errors 2 which can be drizzled as an extensive variable. 
However, caution must be taken when interpreting the resulting 
errors for the following reasons. In principle, if we were deal- 
ing with statistical noise only, by observing the same spot in the 
sky several times, as we effectively do when we combine dif- 
ferent rasters in the drizzle, the error of the combined measure- 
ments would go down. If there is a systematic error it will al- 
ways be the same in sign and magnitude and will not cancel with 
repeated measurements. Furthermore, the errors in the regions 
where the error/signal ratio is better than the calibration accu- 
racy (explained in detailed in Section 3.4) should not be quoted 
as absolute measurement errors. This explains why we have an 
increase of the estimated "errors" where the signal is stronger. 
Therefore, these maps are meaningful only in the "outskirts" of 
the galaxy but not in the starburst region and for this reason we 
masked the starburst regions in the "error" maps shown in the 
right column of Figure [T] In order to give a complete picture of 



the final "noise" levels we reach in each map, we show in Figure 
[TJalso the S/N maps. 



3. Results 

3.1. FIR line intensity and continuum maps 

Figure [2] shows the line intensity maps of the 4 lines in their 
original spatial resolution. All maps are North-South oriented. 
The crosses represent the center of M82. The brightest integrated 
line emissions are very similar at all wavelengths but at the faint 
emission level, some important differences are noticeable. The 
two lines tracing exclusively the neutral atomic gas, namely [OI] 
63 fim and [OI 145] /mi, show a spherical morphology, while the 
[CII] emission, that arises mainly in the neutral atomic medium 
but also from the ionized medium, shows a weak elongation to- 
wards the minor axis. This elongation becomes more evident in 
the [OIII] 88 /mi line intensity maps, that traces purely ionized 
gas. In all maps, emission up to 1 kpc on both sides of the disk 
along the minor axis is detected. From these maps it is evident 
that the ionized medium has an emission elongated along the mi- 
nor axis of the galaxy while the emission arising from the neutral 
medium is more spherically distributed, although still detected 
along the minor axis of the galaxy. 

Figure [3] shows the continuum map derived from the observation 
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Fig. 3. Offline continuum emission obtained from the continuum 
in the [OIII] 88 yum line. Units are in Jy. The black cross corre- 
sponds to the adopted center of M82. 



of the [OIII] 88 yum emission line. The continuum maps at the 
other wavelengths show a similar morphology and therefore they 
are not shown here. The continuum morphology is spherical and, 
at the faint level of the emission, does not show a marked asym- 
metry along the galaxy's minor axis. The brightest continuum 
emission region is displaced ~ 6" West with respect to the cen- 
ter of the galaxy (black cross on the Figure). This is marginally 
larger than the final astrometric uncertainties of our maps (~ 4"). 
The peak emission moves towards the east with increasing wave- 
lengths indicating that the western part is hotter than the eastern 
part in the starburst region of M82. 

3.2. Velocity dispersion maps 

Figure [4] shows the velocity dispersion maps obtained from 
FWHM of the line fits, converted in velocity and corrected for 
the instrumental profile. The overall shape is similar in all lines 
and it is formed by two well distinct kinematic structures: the 
disk and the outflow. In general the two cones of the outflow 
have a velocity dispersion higher (from ~ 50 to ~ 100 km s" 1 or 
a factor from 1.5 to 2) than the disk. 

All four maps clearly show emission in the outflow aligned along 
the minor axis. This enhanced emission originates from a point 
that coincides with the center of the disk (as marked by the black 
cross in the figure). The yellow ellipse shown in the Figure corre- 
sponds to the stellar disk defined by the emission in the ^-band, 
as shown in Figure [9] and explained in Section 3.5. 

3.3. Radial velocity maps 

Figure [5] shows the radial velocity maps, obtained from the 
wavelength corresponding to the fitted lines peak converted into 
velocity, from which we have subtracted the systemic velocity 
of M 82 equal to 203 km s _1 (Shopbell and Bland-Hawthorn 
119981 ). The overall S-shape morphology of the velocity field is 
similar in all lines and confirms previous findings: the east part 
of the disk and the north outflow are receding from us while the 
west part of the disk and the southern outflow are approaching 
us. The region where the velocity is zero is consistent in all 4 
maps and it corresponds to the galaxy's center (black crosses). 
The velocity in the disk goes from — 80 to 100 km s" 1 in all 
lines but the [OI] 145 yum line, which is shifted by — 1-20 km s" 1 
with respect to the others. This shift is present in all raster 



positions suggesting a systematic error (see next section for 
details) rather than a real physical effect. We have also derived 
rotation curves along the major axis of the galaxy (see Figure 
[6]). All curves are consistent with each other (except for the 
velocity shift in the rotation curve obtained with the [OI] 145 
fim line) and they also agree with those shown by Shopbell 
and Bland-Hawthorn (1998 ), especially with those in Ha and 
12 CO(J2 -> 1). 

It is important to note that we do not detect line splitting in 
any of the observed lines. This is particularly intriguing for the 
[OIII] line at 88 yum. Indeed, line splitting is clearly detected 
in Ha in the outflow, above 300 pc from the galaxy plane 
along the minor axis of the ga laxy, with velocity separation 
of ~ 200 km s" 1 (Greve et al. [20041) . well above the PACS 
spectral resolution of the [OIII] line at 88 yum (120 km s" 1 ). It is 
not clear whether line splitting is detected in the optical [OIII] 
lines: Shopbell and Bland-Hawthorn dl998b report marginally 
detected [O III] 5007 A line splitting in some regions of the 
mapped area but the low S/N prevented them to fit double 
Gaussians to their [O III] 5007 A data and therefore to be more 
quantiative on this issue. As will be discussed later, this is not 
the only striking difference we find between the gas emitting in 
the [OIII] line at 88 yum and that emitting in Ha. 

From the [OIII]/[OI] 63 yum map we will present in Section 
3.5, where the bipolar outflow can be best determined, we mea- 
sure an opening angle of ~ 50° for both sides of the outflow, 
very similar to the value found by Walter, WeiB and Scoville 
(2002) in the CO maps. Taking into account the disk inclination 
of M82 (i = 81 ), the outflow opening angle, and assuming a 
simple jet-like geometry, we can derive the mean de-projected 
velocity dividing the observed velocity (~ 40 - 45 km s" 1 ) by 
sin($l° - 50° 12). This results in a mean deprojected velocity 
from ~ 75 to 85 km s" 1 , similar within the wavelength uncer- 
tainties, in both cones of the outflow and at both wavelengths, 
i.e. for both the ionized gas traced by the [OIII] emission and the 
neutral gas traced by the two [OI] lines and the [CII] line. 



3.4. Uncertainties 

The dominant source of flux uncertainties in our data are those 
related to the spectrometer flux calibration, that is -30% in both 
the blue and the red channels (Poglitsch et al. 2010 ). In the raster 
observations presented here, another source of error arises from 
the fact that for all observed lines we used the bright line mode, 
i.e less grating steps per grating scan than in faint line mode. 
The reason for having chosen this observing mode is that these 
maps were executed during the Performance Verification phase, 
to check that the PSF sampling scheme for the mappping mode 
worked as expected. In order to test it we had to observe as many 
wavelengths as possible along the full PACS wavelength range 
(55-220 yum), and faint line mode would have taken a prohibitive 
amount of time. 

In bright line mode each of the 16 spectral pixels of one spatial 
pixel (hereafter spaxel) does not see the full line, as in the case 
of faint line- scan. As a consequence, one of the crucial steps of 
the faint line mode data reduction, bringing the response of the 
16 spectral pixels of each spaxel to the same continuum level 
(the so called FlatField task) cannot be applied in the case of 
bright line mode. Figure [7] shows an example of a spectrum be- 
fore rebinning and combining the nods, for the central spaxel 
of the central raster position at 63 yum. Each spectral pixel of 
the central spaxel has a different color. The figure clearly shows 
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Fig. 4. Velocity dispersion maps obtained from the fitted line widths deconvolved with the instrumental profile width, in km s" 1 . 
The black cross corresponds to the adopted center of M82. The yellow ellipse corresponds to the stellar disk as described by the 
K s -ba.nd emission shown in contours in Figure [9] 



how small is the overlapping spectral range seen by subsequent 
spectral pixels of each spaxel and illustrates the difficulties in flat 
fielding this type of observations. The lack of the flat field could 
result in distorted line profiles, and therefore it could affect the 
velocity peak, the FWHM and the total flux. 
To evaluate the impact of this particular observing mode on the 
intensity maps, we have compared the LWS line flux of M82 
given in Brauher, Dale and Helou (12008b and by Colbert et al. 
(1999) with ours. We have assumed a circular Gaussian LWS 
beam profile, with a FWHM equal to 80" at all wavelengths but 
at 158 jim where we assumed a FWHM equal to 70" (Lloyd 
2003 ). We integrated the line and continuum PACS flux under 
the assumed LWS profiles and we compared the obtained fluxes 
with those observed with LWS corrected for the LWS extended 
source correction factor (see Table 7 of Brauher, Dale and Helou 
2008 ). The results are shown in Figure [8] The LWS and PACS 
fluxes do agree within the respective uncertainties adopted equal 
to 30% for PACS a nd 20% for LWS as claimed in Brauher, Dale 
and Helou (2008 ), with the possible exception of the contin- 
uum flux density at 88 jim which only marginally agrees with 
the adopted uncertainties. We conclude that the bright line mode 
does not introduce uncertainties greater than the nominal PACS 
spectrometer flux uncertainties. 

The nominal uncertainties due to the wavelength calibration 
are ~ 20 km s" 1 and ~ 40 km s" 1 in the blue and red chan- 



nel respectively (Poglitsch et al. 12010b . A distortion of the line 
profile due to the bright line mode could introduce further uncer- 
tainties which is difficult to quantify. However, we expect these 
uncertainties to increase with decreasing line intensity where the 
distortion of the lines profiles become more significant. Among 
the four observed lines, the [OI] 145 fim line is by far the weak- 
est. Therefore, we expect the velocity line peak of this line to 
be more affected than the others due to this observing mode. 
This could in principle explain the shift observed in the rotation 
curve of this line with respect to the rotation curves obtained 
in the other lines and visible in the mono dimensional rotation 
curves shown in Figure [6j On the other hand, the issues with 
the bright-line mode would modify the line profiles in a random 
way, because the misalignment of the spectrum among the spec- 
tral pixels of the same spaxel depends on the transient history 
of each pixel and on the mispointing, which has no systematic 
direction. Since the observed shift in the [OI] 145 jim line is sys- 
tematic in all the raster positions, we believe that the problem 
arises from a wavelength calibration issue at this wavelength. 
We note however, that the fact that the rotation curves in the 
other lines agree with those previously published indicates that 
the bright line mode does not introduce significant additional un- 
certainties even at 145 jjm and that the data are dominated by the 
nominal wavelength calibration uncertainties. 
The deformation of the line- shape due to the bright line mode is 



6 



A. Contursi et al.: Spectroscopic FIR mapping of the disk and galactic wind of M82 with Herschel-PACS. 



[0 I] 63]L/m 



[0 III] 88/vm 



69° 40' 00' 




175.001 



121.67 



68.33 



15.00 



-38.33 



-91.67 



-145.001 



■9° 40' 00' 




09 h 56 m 00 s 



175.00 
121.67 
68.33 
15.00 
-38.33 

-91.67 S 
-145.00 I 



[0 l]145jLvm 



[Cll]158/vm 



o 41' 00" 



69° 40' 00" 




09 h 56 m 00 s 55 m 50 s 



175.001 

121.67 
68.33 | 
15.00 

-38.33 

-91.67 
-145.00 



o 41' 00" 



;9° 40' 00" 




09 h 56 m 00 s 55'" 50 s 

a J2000 



175.00 
121.67 
68.33 
15.00 
-38.33 

-91.67 
-145.00 



Fig. 5. Radial velocity maps obtained from the fitted line profiles at each position. The black cross corresponds to the adopted center 
ofM82. 



also a random process and it would affect the fitted FWHM from 
which we have derived the velocity dispersion maps shown in 
Figure [4] The fact that all maps show the same kinematic mor- 
phology, ensures that the line deformation does not affect the 
overall velocity dispersion in a systematic way, producing fake 
structures. 



3.5. Ratio Maps 

Figure|9]shows the [OI] 63 /mi /[CII], [OIII] /[OI] 63 /mi, [OIII] 
/[CII] and [OI] 145/mi /[OI] 63 /zm ratio maps. Before making 
each ratio map, the two images have been smoothed to the same 
resolution, if necessary. The HST Ha contours (Mutchler et al. 
120071) are overplotted on the [OIII] /[OI] 63 /zm ratio map (see 
discussion in the next Section). 

The ratio maps are very different from each other: the [OI] 63 
/zm/[CII] shows an overall asymmetry along the disk major axis, 
while in the [OIII]/[OI] 63 /zm map and, at a poorer resolution 
also in the [OIII]/[CII] ratio map, the bipolar outflow is very 
clearly delineated. On the [OI] 63 /zm/[CII] ratio map in Figure 
[9] the K^-band contours (Veilleux, Rupke and S waters 2009J de- 
scribing the stellar disk are also plotted. Based on these contours 
we have defined two ellipses: one external (yellow in Figure [9]) 
describing the overall stellar disk and one internal (purple in 
Figure [9]) that defines the brightest part of the stellar disk in- 
cluding the starburst region. 



While the region where the [OI] 63 /mi /[CII] ratio is high 
corresponds well to this internal ellipse, the region where the 
[OIII]/[OI] 63 /zm ratio is higher arises from the north part of 
the starburst region. A comparison with the l2 CO(J 1— >0) map 
(Figure 1 of WeiB, Walter and Scoville l2005l) with the [OIII]/[OI] 
63 /zm ratio image shows that the outflow coincides with the re- 
gions marked as outflow in their figure. The northern outflow 
is a factor of ~ 2 brighter in the [OIII]/[OI] 63 /zm ratio than 
the southern one. In contrast, at optical wavelengths the north- 
ern outflow has not been detected in the [OIII] emission line 
at 5007 A while the southern outflow is clearly visible at this 
wavelength (Figure 2 of Shopbell and Bland-Hawthorn 1998b . 
It seems unlikely that the difference in the observed [OIII] 5007 
A line emission in the two cones of the outflow of M82 is due 
to an electron density higher than the critical density of the 
[OIII] line at 5007 A (6.3xl0 5 cm -3 ) in the northern outflow. 
In fact, the estimated averaged electron density in the outflow 
is < 300 cm -3 (Yoshida, Kawabata and Ohvama 12011b , and the 
average electron density associated with the filaments emitting 
in Ha is ~ 15 cm" 3 (Shopbell and Bland-Hawthorn [T998]) . As 
we will argue in the next section, we believe that this difference 
can be explained with the extinction being higher in the northern 
than in the southern outflow. 

The [OI] 145/[OI] 63 /mi ratio is shown in Figure 9] This ra- 
tio varies only by a factor of 2 all over the mapped area, and 
neither the outflow nor the stellar disk are clearly visible. In 
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the starburst region defined by the purple ellipse, this ratio ap- 
pears slightly weaker than elsewhere. Figure 2 of Tielens and 
Hollenbach (1985 ) shows how this ratio depends on the gas tem- 
perature for various gas densities. For gas temperatures higher 
than -300 K and a fixed gas density this ratio is basically con- 
stant and it decreases with increasing density. The comparison 
with the range of ratios observed in M82 ( ~ 0.05-0. 1 ) suggests 
that the gas temperature is greater than ~ 300 K everywhere and 
that the density is higher in the central region where this ratio 
is lower. We will see in Section 4.2 that more detailed modeling 
confirms that this is the case. 



Finally, Figure 10 shows the line over the continuum as mea- 
sured around the [OIII] and [OI] 63 jim lines. These figures con- 
firm that both the neutral and the ionized gas traced by the FIR 
fine structure line participate in the outflow, but that the ionized 
gas is more collimated than the neutral gas. This fits well in the 
picture where the cold entrained gas forms an outflow less col- 
limated than the hot gas (Ohyama et al. 2007 ). However, the 
ionized gas traced by the [OIII] line at 88 jim is kinematically 
coupled with the neutral atomic gas and it is not kinematically 
coupled with the ionized gas traced by the Ha emission (Greve 
2004) because it moves at a much lower velocity than the hot 
(10 4 K) gas. 



3.6. Comparison with Ha 

In order to understand where the outflow outlined in the 
[OIII]/[OI] 63 fim ratio map is located with respect to the well 
known bipolar outflow of M82 observed at other wavelengths, 
the HST Ha contours are overplotted on the [OIII]/[OI] 63 /im 
ratio map in Figure[9] Overall, the Ha emission is more extended 
than the [OIII]/[OI] 63 /mi ratio (orange contours in the upper 
right panel of Figure [9]). The bulk of the Ha emission, i.e. the 
"eyes" defined in Shopbell and Bland-Hawthorn (1998 ) (black 
contours in the upper right panel of Figure [9]), are displaced with 
respect to the region with of highest [OIII]/[OI] 63 /zm ratio. 
Actually, the region with the brightest [OIII]/[OI] 63 jim ratio al- 
most envelopes the Ha "eyes". Also, in the northern outflow the 
Ha emission is displaced to the west with respect to the region 
with the highest [OIII]/[OI] 63 jim ratio, while in the southern 
outflow both emission are co-spatial. This displacement is most 
likely due to extinction of the optical emission; in fact, the disk 
of M82 is inclined by 81 degree in the sense that its southern 
part is offered to our direct line of sight. Thus, the emission of 
the northern part of the disk and part of the northern outflow 
must suffer from a higher obscuration than that in southern out- 
flow. This agrees very well with what was found by Shopbell and 
Bland-Hawthorn (1998 ) who calculated that, given the inclina- 
tion of M82 disk, the northern outflow should be at least partially 
obscured up to a projected radius of ~ 48". This corresponds to 
the extension of the bright blob visible in the [OIII]/[OI] 63 fim 
map, north-east from the northern Ha emission. 

Extinction could also explain the fact that the northern out- 
flow is detected in [OIII] at 88 /mi but not at 5007 A. We tried 
to make a very rough estimate of the extinction in the follow- 
ing way: we estimated the logarithm of the Ha/Hj3 ratios of both 
sides of the outflow (0.6 and 1.2 for the southern and the north- 
ern respectively) at a distance of ~ 500 pc from the galaxy cen- 
ter from Figure 16 of Heckman, Armus and Miley (1990). With 
these values, we could calculate the A v /R v in the northern and 
southern outflow using the Calzetti law (Calzetti, Kinney and 
Storchi-Bergmann [19941 ). We derived a very crude value of the 
observed [OIII] 5007 A line flux in the southern outflow from 
Figure 2 of Shopbell and Bland-Hawthorn (1998 ) and, using 
a slab geometry, the corresponding dereddened flux. Assuming 
that the intrinsic [OIII] 5007 A emission in the northern out- 
flow is equal to that of the southern outflow, we estimated the 
expected flux, using the A v /R v values obtained for the north- 
ern outflow. We find that the expected [OIII] 5007 A flux is 
~ 2 x 10" 17 erg s" 1 cm -2 for the slab model which is appropri- 
ate for the Calzetti law. This flux is smaller than the minimum 
flux reported by Shopbell and Bland-Hawthorn dl998b in their 
Figure 2, and it is compatible with the non-detection of the [OIII] 
5007 A line in the northern outflow. We stress, however, that this 
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Fig. 8. Comparison between the LWS and PACS fluxes for lines (left panel) and continuum (right panel). The solid line indicates 
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calculation is highly uncertain, first, because of the very uncer- 
tain input values which we had to deduce from the published 
figures, second because there is no physical reason to assume 
that the intrinsic [OIII] 5007 A emission in the two cones of the 
outflow should be the same. 

One important question is whether in the outflow, the bulk of 
[OIII] 88 yum line emitting gas is associated with the bulk of the 
ionized gas emitting emitting in Ha. As already mentioned, these 
two emission components show different kinematics, with the 
[OIII] 88 yum emitting gas being slower (~ 75 km s" 1 ) than the 
Ha emitting gas (~ 600 km s" 1 , Shopbell and Bland-Hawthorn 
1998 ). We have also seen that the [OIII] 88 yum emitting gas does 
not show any line splitting, contrary to what has been observed 
in the outflowing gas emitting in Ha. The Ha line splitting in the 
outflow of M82 has been interpreted as due to the confinement 
of the Ha emitting gas in the walls of the outflow (Greve et al. 
2004). These facts suggest a possible different spatial distribu- 
tion of the ionized gas traced by the Ha and the [OIII] 88 yum 
lines: while the first is confined in the outflow walls, the second 
is more concentrated along the minor axis of the galaxy. 
This may also suggest that these gas components are ionized by 
two different physical processes. In fact, in the commonly ac- 
cepted scenario, the bulk of the Ha emitting gas is confined in 
the walls of the outflow, because it arises, at least in part, from 
shocked gas in the region of interaction between the X-ray emit- 
ting gas and the halo. This naturally explains the observed Ha 
line splitting. Instead, the [OIII] 88 yum line may arise principally 
from ionization due to the starburst photons. Although we do not 
have enough data available to prove that this is the case, previ- 
ous results seem to support this scenario. Shopbell and Bland- 
Hawthorn (119981 ) show that in the southern outflow, the optical 
[NII]/H<2 ratio is lower along the minor axis than further away 
from the axis, probably reflecting a higher ionization. This is in- 
terpreted as radiation produced by the two star forming clusters 
in the disk escaping along the minor axis of the disk. Veilleux, 
Rupke and Swaters ('2009') reached a similar conclusion for both 
parts of the outflow studying the variation of the ratio between 
the near-infrared H2 and the PAH emission. 

We can summarize the main results obtained in this Section 
as follows. Both neutral and ionized atomic gas are detected in 
the outflow and in the disk. In the outflow, the ionized gas emit- 



ting in the [ OIII] 88 jim line is more collimated than the neutral 
gas and also than the ionized gas emitting in Ha, which is more 
confined in the walls of the biconical outflow. The ionized gas 
traced by the [ OIII] 88 \im line and the neutral atomic gas show 
the same kinematics in the disk and in the outflow which differ 
from the kinematic of the Ha emitting gas. These facts suggest 
a scenario where, while the Ha filaments in the outflow arise, 
at least in part, from shocked gas in the regions of interaction 
between the outflowing hot plasma emitting in X-ray and the 
halo of the galaxy, the bulk of the [ OIII] 88 \im line emission 
arises from gas photoionized by the starburst radiation along 
the galaxy minor axis. The comparison between the optical and 
infrared emission suggests that the base of the northern outflow 
is more extincted than the southern one, probably obscured by 
the disk itself. 



4. Analysis 

In this section, we will address the origin of the neutral material 
detected with PACS, with particular focus on its origin in the 
outflow. We have seen that this is certainly material contributing 
to the outflow and the question is whether this emission can 
be explained by classical Photo-Dissociation-Regions (PDRs) 
emission and, if this is the case, how the physical conditions of 
the PDRs in the outflow and in the galaxy disk compare to each 
other. We will also compare the energetics and the kinematics 
of the neutral medium with that of the other gas phases detected 
in the outflow of M82. 



4.1. Photo-dominated regions diagnostics 

The [CII] and [OI] lines are the main coolants of the atomic 
medium. They are collisionally excited by the atomic gas heated 
from electrons ejected from grains via photoelectric effects. 
Since C + has an ionization potential equal to 11.3 eV, it is 
mainly excited by Far UV (FUV) photons (hv< 13.6 eV) es- 
caping from HII regions (but it can also be excited in the HII 
regions themselves as we will discuss in Section 4.1.3). The 
regions where FUV photons dominate the ISM physics are 
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Fig. 9. [01] 63 yum/[CII], [OIII] 88 /mi /[OI] 63 /mi, [OIII] 88 yum/[CII] and [01] 145 /rni/[OI] 63 /mi ratio maps. The K,-band 
contours are overplotted on the [01] 63 yum/[CII] ratio map. They were used to define the outer ellipse tracing the stellar disk and 
the inner ellipse that delineates the starburst region. Ha contours are overplotted on the [OIII] 88 yum /[OI] 63 yum ratio map. The 
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Fig. 10. [OIII] 88 yum /continuum (left) and [OI] 63 yum/continuum (right) ratio maps. The continuum is from the baseline level of 
the line fit. 



called Photondominated or Photodissociation Regions (PDRs, 
Hollenbach and Tielens 1999 and reference therein). Since the 
atomic gas cooling is mainly traced by the [CII] and [OI] line 
at 63 yum and the heating is due to the photoelectric effect on 
the grains, the ratio between the total emission from the FIR 
fine structure lines and the total infrared (TIR) emission from 



grains gives the total photoelectric yield. Comparing the ob- 
served yield with that predicted by the model one can derive 
fundamental physical parameters of the neutral ISM associated 
with the PDRs. 

In this section, we will show the results obtained by apply- 
ing the PDR model from Kaufman et al. (fT999b to the PACS 
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data of M82. By applying the model to each pixel of our maps, 
we are able to produce, for the first time for an extragalactic 
object, high resolution, fully sampled maps of the FUV interstel- 
lar radiation field (ISRF) expressed in Go unit^ of the neutral 
gas density nu, of the gas temperature at the surface of the mod- 
eled clouds T s and of the cloud beam filling factors (p. Because 
of different step sizes in mapping among all wavelengths and 
final dimensions of the various continuum maps, we rebinned 
all maps to the grid used for the longer wavelengths mapping 
(6 // /pix), reduced all images to the same size, smoothed all maps 
to the coarser spatial resolution of the [CII] map, then aligned 
all to the [CII] emission peak. In the following analysis we will 
assume that the [CII] and [01] lines arise from PDRs also in the 
outflow, and we will check a posteriori if this makes physically 
sense. However, in order to correctly and successfully apply the 
PDR diagnostics one has to make some corrections to the ob- 
served fluxes that we will explain in detail in the following sec- 
tions, before showing the results. 

4.1 .1 . Total infrared emission map calculation 

One fundamental input to the PDR modeling is the total IR emis- 
sion (TIR from ~ 30 to 1000 yum). Since our aim is to derive a 
solution for every pixel of the M82 mapped area, we need to have 
also the TIR emission in each pixel. But the IR continuum data 
available from previous observations (IRAS, ISO, SPITZER) do 
not have the spatial resolution that matches that of PACS at these 
wavelengths. Fortunately, thanks to the brightness of M82, we 
could use the continuum detected aside the targeted lines and in 
the corresponding parallel channel, as photometric points from 
which we could derive the FIR ( from ~ 40 to 200 yum) emis- 
sion at each pixel. This gives a total of 8 points per pixel at 63, 
72.8, 78.9, 145, 157, 176.8 and 189.6 yum. We fit the continuum 
fluxes with a modified blackbody, in the optically thin emission 
approximation and with a fixed value of the grain emissivity, i.e. 
with ft = 1 in each pixel. However, the agreement of the set of 
continuum data points obtained from the blue band lines and the 
corresponding red parallel points do not agree very well with 
those obtained from the continuum aside the red lines and the 
corresponding blue parallel data. In particular, the data derived 
from the parallel channels of the [CII] and [OI] line at 145 yum 
(i.e. 72.8 and 78.9 yum) are often systematically higher than the 
data values at the other wavelengths. The reason for this is not 
well understood, but it might correspond to an higher continuum 
distortion at these wavelengths introduced from the bright line 
mode with respect to other wavelengths which increases the un- 
certainties. 

Therefore we restrict the fitting to the 4 photometric points cal- 
culated on the blue lines and the corresponding parallel data. In 
this way we include in the gray body fitting the point at the high- 
est possible wavelengths, allowing a better determination of the 
far infrared emission. An example of such a fit is shown in Figure 
[TTJfor one pixel. From the black body fit, for each pixel we cal- 
culated the flux density at 60 and 100 yum which is necessary for 
the comparison with diagnostics used in the pre-Herschel era, 
and from these flux densities we derived the FIR (from 42.5 to 
122.5 jim) emission using the Helou (1988) formula : 

FIRiWIm 2 ) = 1.26 x 10" 14 x (2.58 x I 6 ^ m (Jy) + /ioo^(/y)).(l) 

The blackbody fitting also results in a dust temperature for 
each pixel that we used to build the map shown in Figure [12 
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Fig. 11. Example of a modified blackbody fit on one pixel in the 
center of M82. The points correspond to the continuum averaged 
detected with PACS on both sides of the [OI] line at 63 yum and 
[OIII] at 88 yum and on the corresponding parallel channel. 



where also the Ha contours and the smaller of the two ellipses 
shown in previous Figures are shown. Taking into account the 
flux uncertainties we estimated an uncertainty on the dust tem- 
perature of + 3 K. The dust temperature image shows very 
clearly the bright central region included in the ellipse and cor- 
responding to the brightest Ha emission. Here the dust temper- 
ature ranges from ~ 45 to ~ 55 K, typical for a starburst re- 
gion. In the map is also visible an extention of the dust emission 
above the disk along the minor axis which coincides with the 
Ha extension along this direction and with the outflow seen in 
the [OIII]/[OI] 63 yum ratio map. The dust in the southern out- 
flow is slightly hotter than in the northern outflow. 
From the FIR emission we calculate the TIR emission needed as 
input to the Kaufman et al. (1999 ) model by using eq. 3 given in 
Dale et al. ( 2001 ) which calculates the conversion from FIR to 
TIR from the 60/100 /mi ratio. 

4.1 .2. Geometrical effect correction factor. 

The PDR model we consider in this work assumes the emission 
is generated by a one-dimensional semi-infinite slab illuminated 
from one side. The model also considers [CII], FIR and [OI] at 
145 //m optically thin, and [OI] at 63 fim optically thick, such 
that in a model where a cloud is illuminated from one side, only 
one side of the cloud emits in the [OI] line at 63 yum while [CII], 
FIR and [OI] at 145 yum emits from both the front and the back 
of the cloud. At the PACS resolution each beam includes many 
clouds, and the orientation of the illuminated side will be ran- 
domly distributed. This is equivalent to consider that the beam 
includes one cloud illuminated from all sides. The model would 
then predict half of the [OI] at 63 jjm. This is why we multi- 
ply the observed [OI] at 63 yum by 2 before applying the PDR 
models. Since the [OI] at 145 fim is considered to be thin in the 
model, we do not have to apply any correction to it. 

4.1 .3. Contribution to the [CII] emission from the ionized gas 

[CII] can also arise from the ionized medium. Since the PDR 
model considers [CII] arising only from the neutral medium, be- 
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Fig. 12. Dust temperature map. Units are °K. Ha contours are 
overplotted. Note that these contours are at a resolution poorer 
than what is shown in Figure [9] to match the pixel size of the dust 
temperature map. 



fore applying the model to the observations we have to subtract 
this contribution from the observed [CII] emission. The percent- 
age of the [CII] emission arising from ionized gas with respect to 
the total observed emission, is minimum ( ~ 10%) in dense HII 
regions, since here most of the carbon is double ionized, but it in- 
creases significantly in the diffuse ionized medium, the precise 
amount depending on the physical parameters of the medium 
(Nagao et al. 1201 lb . In external galaxies the beam typically en- 
compasses regions with very different physical conditions, there- 
fore one has to guess the average properties of the observed 
medium in each spatially resolved region, in order to estimate 
the [CII] contribution from ionized gas. This will obviously in- 
troduce some uncertainty in the derived parameters. 

The most straightforward way to determine the fraction of 
[CII] arising from the ionized medium is to use the scaling re- 
lation between the ( 3 P i-> 3 P ) [Nil] line at 205.178 jim and the 
[CII] line (Abel 2006 ). Since these have similar critical electron 
densities, (80 and 50 cm -3 respectively), their ratio is almost in- 
dependent on the electron density. Unfortunately, we cannot use 
this method because we did not observe M82 in the [Nil] 205 
fim line with PACS. A less direct, and more uncertain, way to 
calculate the contribution from ionized gas to the observed [CII] 
line flux, is to use the [Nil] 122 fim line. Since [Nil] 122 /mi 

has a higher critical electron density (3.1xl0 2 cm" 3 ) than that 4 - 2 - Results from the PDR modeling 



VCU^Xed = VCIUZse HII = LI >< WIT* . (2) 

In a recent work, Croxall et al. (2012), show how the 
[CII]/[NII] 122 ratio varies with electron density (their Figure 
11). We do not have a direct electron density measure in the 
outflow but we know from previous works that the average 
electron density of the ionized gas traced by the Ha is ~ 15 cm -3 
(Shopbell and Bland-Hawthorn 1998) and the estimated aver- 
aged electron density in the outflow is < 300 cm -3 (Yoshida, 
Kawabata and Qhvama [20111) . Since we do not know a priori 
from which of these these media the [CII] arises, we have 
calculated the amount of [CII] arising from ionized gas for 
these two cases by using the ratio published by Croxall et al. 
(I20T21) . From this work we derive a [CII]/[NII] 122 ratio equal 
to 3.7 and 0.7 for n e = 15 and n e = 300 cm -3 . Note that these 
two density values correspond more or less the minimum and 
the maximum correction possible, in order not to have negative 
[CII] results. In each case, the resulting [CII] fluxes were then 
subtracted from the observed [CII] map, in the region defined 
as outflow. The resulting percentages of the [CII] emission 
arising from the ionized gas are ~ 10, 30 % in the disk and in 
the starburst, and ~ 11 (47) % and ~ 11 (40) % in the north 
and south cones of the outflow based on n e = 300 (15) cm" 3 
respectively. 

We stress, however, that these estimates of the contribution 
to the observed [CII] flux from pure ionized medium are very 
uncertain especially when applied to regions with an unknown 
and complex mix of different ISM phases. Moreover, because of 
the lack of a fully sampled map in the [Nil] 122 jim line map, 
we were forced to use an average value of the [Nil] 122 emis- 
sion for each of the four macro-regions we have defined in the 
galaxy. It is for this reason that we decided to run the model also 
with and without the corrections to the observed [CII] emission 
line we have explained above, to check whether the main results 
we derive from the modeling are independent of the correction 
applied. 

Moreover, in principle we could have used as input for the PDR 
modeling the two [OI] lines, [CII] and TIR. We have tried this 
but the uncertainties in each of the line is such that the result is 
much nosier than using 3 input values with only one [OI] line at 
the time. Therefore, we run the model using always only three 
input values, i.e the [CII], the FIR and one of the [OI] lines at the 
time, and with and without the ionized contribution correction to 
the [CII] emission line (i.e. 6 input sets in total). 



of the [CII] line, the ratio of these lines depends strongly on the 
density. 

We have observed M82 in the [Nil] 122 /mi line with three 
PACS spectrometer pointings centered on the disk, and on the 
north and south outflow as shown in Figure 13 We have divided 
the galaxy in four regions: the starburst, the north and the south 
outflow and the more diffuse disk and we have calculated the 
mean [Nil] emission for each of them using the average of the 
spectra belonging to each of the defined regions. The resulting 
average [Nil] intensities in the starburst, disk, and in the northern 
and southern parts of the outflow are 87.4, 1.0, 5.3 and 3.4 x 
10" 17 W m" 2 . We assume that in the disk and in the starburst 
region of M82 the [CII] fraction arising from ionized gas scales 
with the observed [Nil] following the relation valid for the dense 



By applying the PDR model for each input data set, we obtained 
Go and n H for each pixel where all three input values are de- 
tected. From Go and n H we then derived the corresponding gas 
temperature at the surface of the modeled cloud (T s ) and the area 
beam filling factor (0), i.e. the percentage of the beam covered 
by the emitting clouds, obtained by dividing the [CII] emission 
corresponding to the Go and n H solutions and the input [CII] 
value. Since this can be done for each pixel we were able to ob- 
tain 6 such maps for each input dataset defined as explained in 
the previous section. We stress, however, that given the flux un- 
certainties and the fact that a model is an oversimplification of 
much more complex real physical conditions, we estimate the 
uncertainties of the derived PDR model parameters to be a fac- 
tor ~ 2. Also we do not really trust the pixel-to-pixel variation 
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Table 1. Main PDR physical parameters obtained by modeling 
the [CII], [01] 145 and FIR emission with and without correct- 
ing the observed [CII] line flux from the contribution from vari- 
ous ionized gas (see Section 4.1.3 for details). The listed values 
are calculated by averaging the maps shown in figures from 14 
to [T7] over the four macro regions we have divided the galaxy 



in: starburst, disk, north and south outflow (see figure [18] for a 
schematic visualization of these regions). Last column lists the 
averaged value for the [CII] opacity maps obtained as explained 
in Section 4.3. 
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in the PDR solution maps, simply because the model uses a dis- 
crete grid, and given the flux uncertainties it is possible that a 
resulting solution is higher or lower than the solution it would 
have been assigned with smaller uncertainties. Sometimes this is 
visible in a hot/cold pixel in a middle of a well defined morpho- 
logical structure, and we corrected for these interpolating with 
the neighborhood pixels. 

Figures 14 to [17] show the 6 maps obtained using the six sets 
of inputs, for each of the above mentioned physical parameters. 
In detail: the top panels show the solutions obtained using as in- 
put to the PDR modeling [OI] 63, [CII] and FIR. Two of these 
maps were obtained by correcting the observed [CII] emission 
from the contribution arising from ionized gas assuming electron 
densities equal to 300 and 15 cm -3 in the outflow and applying 
equation 2 in the disk and starburst regions. The third solution 
was obtained without correcting the observed [CII] . The bottom 
panels of figures [14] tc fTT] show the three sets of solutions ob- 
tained in a similar way but using as input to the model the [OI] 
line at 145 yum instead of the [OI] line at 63 yum. 

All maps are very similar in morphology except for the gas 
density maps (Figure[l5]). In particular, this figure shows that the 
gas density is enhanced in the starburst region for all the solu- 
tions obtained using the [OI] line at 63 yum, and in the northern 
part of the starburst when the [OI] line at 145 yum is used. 
We find that this difference is most likely due to misalignment 
between the [OI] 63 yum and [CII] images. We explain this in 
details in the Appendix. Once this misalignment has been cor- 
rected, the morphologies of the images obtained by the PDR 
model ing fo r all parameters, included the density n H , agree (see 
, and the values agree within a factor of - 



the fact that the [OI] 63 fim optical depth can significantly vary in 
the galaxy, we will use for the following analysis only the results 
obtained by using as input the [OI] 145 yum line. This reduces the 
uncertainties of our analysis to those related to the [CII] correc- 
tion for ionized gas, since the [OI] 145 fim line is optically thin. 
Moreover, the pointing variation and the the processing affects 
the [CII] and [OI] 145 yum maps in a similar way, since they 
have been observed in the same AOR. The main uncertainty left, 
i.e. that due to the [CII] correction for the ionized gas contribu- 
tion, can be assessed by examining the results obtained from the 
corrected and the not corrected [CII] values. 

Therefore, we now concentrate on the results shown in the 
bottom panels of figures [T4|to[T7| The first thing to notice is that, 



for each of the derived maps, there is no significant difference in 
the morphology among the solutions obtained with no [CII] cor- 
rection, and those with the two different corrections we have ap- 
plied to the observed [CII] flux. We have calculated the average 
values of Go,n H , T gas and cp in the four regions we have consid- 
ered for the different correction to the observed [CII], for each 
case, and listed them in Table[T] The ranges of values of the PDR 
solution parameters are all comparable within the uncertainties 
(a factor of 2~3). We are thus confident that the correction we 
have applied to the observed [CII] does not really play a major 
role. 

As already mentioned, the gas density maps, show a density en- 
hancement in the upper part of the starburst region, correspond- 
ing also to the brightest emission in the [OIII]/[CII] ratio map 
shown in Figure [9] 

The gas temperature is everywhere greater than ~ 300 K as we 
had already deduced from the constancy of the two [OI] lines 
ratio map. The brightest emission outlines the disk of the galaxy 
and no structures corresponding to the outflow are visible. 
The PDR area filling factor maps are very similar in all cases 
and show that there is an enhancement in the galaxy disk with 
a slight asymmetry toward its northwestern part. The area fill- 
ing factors in the outflow are significantly smaller than in the 
disk and very similar to the diffuse more spherically distributed 
emission visible in all line intensity maps. 
One interesting result common to all solutions is that we mapped 
an area sufficiently large to detect relatively diffuse emission 
from the galaxy disk, with Go ~ 100 and nu ~ 50 cm -3 . 
Previous space missions have not had the sensitivity or spatial 
resolution to detect FIR fine structure lines from the low density 
cold atomic medium of galaxies. As demonstrated by our results, 
Herschel will greatly advance the study of the multi-phase ISM 
in external galaxies in the FIR. 



Figure |AJJ ), and the values agree within a factor of ~ two for all 
maps, except the density map for which the biggest difference is 
close to a factor of ~ 3. 

Because of the problems and uncertainties discussed above, and 



Figure 18 shows the logarithmic plot of Go versus n H for 
all solutions obtained with the input sets which included the 
[OI] line at 145 fim. Similar plots derived using the results 
obtained with the re-aligned [OI] 63 fim map, are consistent 
with those shown in Figure 18 The macro regions defined for 
the different corrections to the [CII] total flux are sketched in the 
top right panel. The same color code is used in the three plots 
shown in the figure. The upper panel shows the results obtained 
using as input to the model the observed [CII] emission without 
any correction; the bottom left (right) panel shows the results 
obtained assuming that the electron density of the ionized gas 
producing [CII] to be subtracted to the observed flux is is equal 
to 300 (15) cm -3 . The similarity between all relationships is 
evident, reinforcing our assertion that the correction to [CII] 
does not significantly impact our results. This figure shows 
that n H and Go are correlated and that the correlation agrees 
with what was previously found by Malhotra et al. (120011) by 
modeling a sample of local star-forming galaxies observed 
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Fig. 14. The PDR solution maps obtained for the Far UV Interstellar Radiation Field in units of Go. Left to right are the results 
obtained correcting the [CII] flux from the ionized contribution obtained by considering the densities of the ionized gas equal to 300 
(left) and 15 cm -3 (middle) and with no correction (right). The top raw shows the results obtained using as input [01] 63, [CII] and 
FIR; the bottom raw shows the results obtained using as input [01] 145, [CII] and FIR (see Section 4.3.1 for details). 
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Fig. 13. The PACS spectrometer footprints of the 3 pointings ex- 
ecuted in [Nil] 122 yum are overplotted on the [CII] intensity map 
of M82 (the same as shown in Figure[2]). Note that the size of the 
circles does not represent the size of a spaxel on the sky. 



with the Infrared Space Observatory (ISO-LWS) and that are 
shown as yellow asterisks in Figure 18 Also shown in this 
Figure is the value obtained from Colbert et al. dl999b from 
ISO-LWS observations of M82. We find that the solutions for 
M82's starburst region corresponds well to the PDR solution 



determined by Malhotra et al. and Colbert et al. from the ISO 
data. The ISO-LWS observations are single pointings centered 
on the galaxy centers (i.e do not encompass the entirety of the 
disk) and therefore are biassed towards the warm and dense 
central regions. This explains why these measurements occupy 
the high Go~n H part of the diagram. The results of this work 
extend the correlation by almost 2 orders of magnitude to much 
lower Go-nn. All four components of M82 that we considered, 
follow the same correlation in Go - nu as normal galaxies 
suggesting a common origin of the neutral gas emission for 
all regions, including the outflow. This means that the neutral 
gas as traced by its main coolants is consistent with arising by 
classical PDRs also in the outflow. 



4.3. Column density and [CII] optical depth of the atomic 
neutral medium associated with the PDRs 

The averaged cloud column density of the hydrogen nuclei asso- 
ciated with the [CII] emission, N[cn](H), can be calculated using 
the following equation (Crawford et al. J1985b : 



N [cm (H) 



4.25 x 10 



,20 



1 + 2 x e^~ 92 ^ + (n crit /n H ) 



2 e (-92/T) 



jPDR 
'[CII] 



(3) 



where N[cn](H) is in cm -2 , m erg s 
is the area filling factor, i.e. the fraction of the beam filled with 
clouds emitting [CII],^(C) is the [C + ]/[H] gas -phase abundance 
ratio that is equal to 3xl0~ 4 (Crawford et al[T985 ), n cr i t is the 
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16. The PDR solution maps obtained for the temperature of the gas T k 
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Fig. 17. The PDR solution maps obtained for the area beam filling factor ( 
clouds. 



i.e the percentage of the beam covered by the emitting 



critical density for collision of C + with H atoms equal to 4x 10 3 
cm -3 , T is the gas temperature in K and n H is the gas density. 
Using the observed (and corrected) [CII] emission, the tempera- 
ture, the gas density and filling factor values resulting from the 
PDR modeling we are able for the first time to derive the aver- 
aged cloud column density map of the HI gas associated with the 
PDRs emitting the observed [ CII] in M82. 
The resulting column density ranges from ~ 1 x 10 21 cm -2 up 
to ~ 1 x 10 22 cm -2 for all input data sets except for the ones 
obtained using the [OI] line at 145 jim and [CII] not corrected 
for the ionized gas contribution, for which the column density 
reaches few 10 22 cm -2 . We also derived the [CII] optical depth 
map by applying the following formula (Crawford et al. 1985 ): 



T [CII] 
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The results obtained are shown in Figure [19] The morphol- 
ogy of both the average cloud column density (not shown here) 
and the [CII] opacity maps, are obviously similar. All maps 
have [CII] optical depths much less than 1. The last column of 
Table [T] lists the average T\cii] values of the starburst, the disk , 
the north and south part of the outflow. The map obtained with 
no correction to the observed [CII] line flux shows an opacity 
higher (up to ~ 0.1) than those obtained correcting the [CII] 
observed flux from ionized gas contribution, as expected, while 
the morphology is very similar in all three cases. Figure 19 
shows that the opacity is higher along the disk of the galaxy 



with a decrease toward the starburst region, where the brightest 
Ha emission is also located. This opacity decrease could be 
the result of material cleared out from the winds powering 
the outflow. The opacity drop is even more pronounced in the 
outflow regions, especially in the northern one. 

Some estimates of the column density associated with the 
[CII] line emitting gas already exist for M82. In particular we 
compare our results with those obtained from Crawford et al. 
(1985). These authors derived t\cii\ £ 0.2 and column density 
Ncn > 9 x 10 17 cm -2 by using equation 3 and 4 in the limit 
of high density (n » n cr i t ) and high temperature ( T » 92 K). 
We have calculated a corresponding column density and optical 
depth averaging the values obtained from our maps, in a region 
of M82 which corresponds, as close as possible, to the position 
and beam size used in the Crawford observations. We find an 
average [C + ] column density a few 10 18 cm -2 and a opacity ~ 
0.03, i n good agreement with the limits given by Crawford et al. 
( 1985 ) on these parameters. 



4.4. Outflow energetics 

We can estimate the minimum mass of the hydrogenic nuclei as- 
sociated with the [CII] emitting gas using equation 1 of Hailey- 
Dunsheath et al. (l20T0l) . 

M[cii](H) = x ( OTL^ x (1+2*™+^ 9 



2 e (-92/T) 



In what follows, we calculate the mass, the kinetic energy 
and the mass outflow rate associated with the neutral atomic gas 
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Fig. 18. PDR solutions for M82 obtained by using [CI], [01] at 145 fim and TIR as input parameters to the model. On the right of 
the top panel a schematic view of the masks used to define the starburst (green), the disk (black), and the north (cyan) and south 
(red) outflow is shown. The same color code has been used for the plotting. Also overplotted for comparison, the solutions obtained 
for a local sample of normal star forming galaxies from Malhotra et al. 2001 (yellow stars) and with the ISO-LWS single pointing 
flux of M82 (Colbert et al. 1999). Top: PDR solutions obtained by using as input the [CII] emission not corrected for the ionized 
gas contribution. Bottom left : the results obtained by using the [CII] corrected for the ionized gas contribution assuming a electron 
density equal to 300 cm -3 . Bottom right: the results obtained by using the [CII] corrected for the ionized gas contribution assuming 
a electron density equal to 15 cm -3 . (See section 4.1.3 for details). 



emission obtained in this work (Tabl^2]) and we compare them 
to the corresponding values of the cold molecular gas traced 
by the CO emission, the warm molecular gas traced by the H 2 
emission and the ionized gas emitting in Ha (Table [3]). 
Using equation 5 we calculate the neutral atomic gas masses 
associated with the outflow for each input data set involving the 
[OI] line at 145 jim only (listed in Column 1 of Table [2]). The 
overall uncertainties for these quantities are comparable with 
those quoted for the PDR parameters. 

The total atomic mass in form of PDRs in the outflow is 
2 - 8 x 10 7 M Q depending on the input data, a factor ~ 4- 
15 times lower than the cold molecular mass entrained in the 
outflow as measured by Walter, WeiB and Scoville ( 120021 
M ~ 3.3 x 10 8 Mq) and 3 orders of magnitude higher than the 



warm molecular gas (Veilleux, Rupke and S waters 120091) . Our 
values should be regarded as lower limits because the area on 
which we have calculated the outflow masses, is not defined by 
the intensity of the emission but by the extent of the mapped 
area. In other words, there may be still material further out that 
we have just not mapped. Moreover, this is not the total atomic 
mass carried in the outflow but only that associated with the [CII] 
emitting gas. Therefore the atomic and molecular gas in the out- 
flow could be considered as almost comparable. For comparison, 
the mass of the warm filaments detected at optical wavelengths 
is 5.8 x 10 6 M (Shopbell and B land-Hawthorn [T555J. 

The kinetic energy associated with the neutral gas is Ejy n = 
M[cir\(H) x v 2 where v 2 is the deprojected mean velocity of the 
outflow calculated as explained in section 3.3. We obtained val- 
ues in the range of ~ 1 - 5 x 10 54 erg (Table|2) Column 2) for the 
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two cones of the outflow. From Table [3] we see that these values 
are 1 order of magnitude lower than what is found from the out- 
flowing molecular gas (Walter, WeiB and Scoville l20021) and Ha 
gas (Shopbell and Bland-Hawthorn l 1 9981) equal to ~3xl0 55 erg. 
These values are also ~ 2 orders of magnitude higher than the ki- 
netic energy of the warm H 2 molecular gas calculated by assum- 
ing that the velocity of the warm molecular gas is equal to that 
of the CO molecular gas (Veilleux , Rupke and S waters [20091) . 
Therefore, although not energetically dominant, the atomic gas 
in the outflow still represents a significant contribution to its 
overall energetic. 

We define the dynamical time of the outflow as the time nec- 
essary for its material to travel at a mean velocity equal to Vdeproj 
up to a distance of 450 pc away from the disk along the minor 
axis. We use this length because this is the extent of the mapped 
area above and below the galaxy's disk. We can apply this sim- 
ple calculation to both the neutral and the ionized gas, the latter 
being traced by the [OIII] line emission. The dynamical time for 
both cones of the outflow is ~ 5 x 10 6 yr for the neutral and the 
ionized gas, basically because they all have comparable depro- 
jected mean velocities in the outflow. This time is also compa- 
rable to that estimated for the warm molecular H 2 by Veilleux 
, Rupke and Swaters (120091) for the same distance we have as- 
sumed in our calculation. 

Finally we can calculate the mass outflow rate of the neutral 
gas in the outflow dividing the masses by the dynamical times. 
For each input data set, the mass outflow rates (Column 3 of 
Table [2]) are similar in the northern and southern cones of out- 
flow. The range of the total neutral gas mass rate is 10-25 M© /yr. 
However, if we add the mass outflow rate of the molecular com- 
ponent in the outflow (33 M©/yr, Walter, WeiB and Scoville 
2002 ) the total PDR outflowing mass rate in the outflow is com- 
parable to the total Star Formation Rate (SFR) of M82 (~ 25 
M@/yr ) estimated by Forster-Schreiber et al. (2003 ) giving a 
mass load ( defined as Moutflow/SFR) ~ 2 similar to that mea- 
sured in another starburst galaxy NGC253 (Sturm et al. 1201 It . 



4.4.1 . Comparison amongst the various phases of outflowing 
material 

We can summarize the ISM phases participating to the M82's 
outflow as follows: the hot (10 7 K) gas emitting in X-rays; the 
ionized cooler gas (10 4 K) emitting in Ha; the cold and warm 
molecular gas (Walter, WeiB and Scoville [2002, Veilleux , 
Rupke, and Swates I20091 the cold ( ~ 300 - 400 K) neutral 
gas (this work) and the dust (Alton, Davies and Bianchi 1999, 
Leeuw and Robson 2009 , Kaneda et al. 120101 Roussel et al. 
2010). These components are not all dynamically coupled to 
each other. The cold molecular gas traced by the CO emission 
and the neutral gas traced by the FIR fine structure line emission 
studied in this work have similar averaged velocities and open- 
ing angles which suggests they are coupled. These components 
are decoupled from the ionized gas emitting in Ha which has 
much higher velocity (Greve 2004 ). If we assume that also the 
PAHs emission observed by Engelbracht et al. (12006b and the 
warm molecular gas observed by Veilleux , Rupke and Swaters 
(2009), for which we do not have kinematic information, are 
also coupled with the [CII] emitting gas in the outflow, we have 
gathered all fundamental components of the classical PDRs. 
Probably, the cold gas material is entrained from the disk into 
the outflow by the fast hot wind, in the form of clouds. Once in 
the outflow, the UV photons in the wind itself and/or directly 
from the starburst region, excite the external shells to these 
clouds which emit in the observed FIR fine structure lines, the 



aromatic bands from PAHs and warm molecular gas tra ced by 
H2 in the near infrared. Micelotta, Jones and Tielens (2005 ) 
have shown that PAHs can survive the destruction by collisions 
with thermal electrons in the hot winds only if they are in 
form of cold cloudlets. This supports the idea of having the 
molecular and cold neutral medium in the outflow organized 
in small clouds dragged by the wind, dynamically decoupled 
from the ionized gas traced by Ha emission, but almost equally 
energetically important. 

The analysis of the behavior of the derived Go along the 
minor axis of M82 can help test this hypothesis. In Figure 20 
we show the values of Go we have obtained from the PDR 
modeling along the minor axis of the galaxy as a function of the 
distance from the galaxy's center, for the different input data 
sets we have used. The two panels show the cut for the southern 
(left) and the northern (right) outflow. 

We compare these values to the ISRF values provided by the 
starburst assuming pure geometrical dilution without any extinc- 
tion. This corresponds to the maximum possible contribution 
of the starburst light in the outflow. In order to calculate it we 
assume the following: 1) the starburst has a diameter of 500 
pc (Forster Schreiber et al. 12003b , which corresponds to -30 
", and therefore we consider only the points at distances >15", 
along the minor axis on each side; 2) the ISRF at this distance 
in both cones of the outflow, is equal to what we have calculated 
from the PDR modeling at the same distance. The resulting 
values are indicated with the blue solid line in both panels 
of figure [20] This figure shows that the ISRF values derived 
from the PDR modeling agree within the uncertainties with the 
blue line, except for the Go values obtained by using as input 
to the PDR models the [CII] corrected from the ionized gas 
contribution for a density equal to 15 cm -3 , which is a rather 
extreme case. Agreement with the maximum Go obtained by 
neglecting extinction would indicate that most of the starburst 
FUV photons can travel the outflow, with a low probability of 
being absorbed. Since we know that there is dust in the outflow, 
the only way this can happen is that there are big dust-free 
holes or, in other words, that most of the existing dust must be 
concentrated in clouds with low filling factors. This supports 
the cloudlets interpretation we have proposed above and it also 
shows that the starburst light is sufficient to heat most of these 
cloudlets and that no significant extra heating (i.e. in situ star 
formation) is required. 

We have already discussed in Section 3.6 that the ionized gas 
traced by the [OIII] line at 88 jim and that traced by the Ha emis- 
sion are not kinematically coupled and probably not even co- 
spatial in the outflow. We have also shown evidences that these 
two lines likely trace gas ionized by two different physical pro- 
cesses: the Ha emission arises from a combination of photoion- 
ization by the starburst and shocks: the first process dominates 
the emission towards the outflow axis, the second process traces 
predominantly the gas shocked when the hot flowing plasma en- 
counters the halo, and it is therefore more confined in the outflow 
walls. The [OIII] 88 jim emission is more concentrated along the 
minor axis of the galaxy and it is most likely dominated by pho- 
toionization. Taking also into account that the [OIII] 88 fim emit- 
ting gas flows at the same velocity of the neutral atomic gas, we 
propose that this is gas surrounding the PDR-molecular cloudlets 
ionized by the starburst light. This ionized component though is 
more collimated than the neutral gas component, suggesting that 
the ionized gas associated with the PDR clouds, is decreasing 
towards the periphery of the outflowing material. This is prob- 
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n e =300 [O I] 145jL7m n e =15 [O I] 145jL7m No correction [O I] 145jL7m 




Fig. 19. The [CII] opacity maps obtained by using equation 4. Only the result of the modeling obtained using as input data set [CII] 
(with and without corrections), FIR and the [01] line at 145 /mi are shown. 

Table 2. Atomic gas characteristics in M82 outflow 



input set 


M [cm (HI) 




M 




10 7 M 


10 54 erg 


M /yr 


[01] 145 /im, [CII] corrected n e = 15 cm" 3 








North 


2.1 


1.5 


4.0 


South 


2.9 


1.7 


5.1 


[01] 145 fjm, [CII] corrected n e = 300 cm" 3 








North 


4.3 


3.2 


8.4 


South 


8.4 


5.0 


14.8 


[01] 145 yum,[CII] NOT corrected 








North 


4.1 


3.1 


8.1 


South 


6.1 


3.6 


10.7 



Table 3. Masses, kinetic energies and mass outflow rates in the outflow of M82 for the cold, warm and ionized gas traced by the Ha 
emission, obtained in other works. These values should be compared with the values obtained for the neutral atomic emitting gas 
calculated in this work and listed in Table 1 . 



Gas phases 


M [cm (HI) 


Ekin 


M 


Reference 




M 


erg 


M Q /yr 




Cold Molecular gas 


3.3xl0 8 


3.3xl0 55 


33 


Walter, WeiB and Scoville 2002 


Warm Molecular gas 


1.2xl0 4 


10 51 


0.001 


Veilleux, Rupke and Swaters 2009 


Ionized gas traced by Ha 


5.8xlO b 


2.0x10^ 


3.6 


Shopbell and Bland-Hawthorn 1998 


Neutral atomic gas 


2-8xl0 7 


l-5xl0 54 


10-25 


This work 



1000 r 



100 r 



South 




(J) SB geometrical dilution 

+ [CII] correction n = 15 cm 
H [CII] correction n = 300 cm 
% No [CII] correction 



1000: 7 



100; r 



North 




10 



15 



23 



35 



53 



8D 



15 



23 



35 



53 



80 



Distance perpendicular to the disk (") 



Fig. 20. The observed FUV ISRF in Habig units along the minor axis of M82 are plotted as a function of the distance from the 
galaxy center in ", for the south (left) and north (right) outflow separately. The blue symbols connected by the blue line, are the 
values calculated in the pure geometrical dilution hypothesis. The starburst region, which occupies the central 250 pc (~ 15 ") is 
not included. 
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ably due to the fact that the number of UV photons necessary 
to excite the [OIII] line (> 35.12 eV) directly emitted from the 
starburst decreases faster than the number of the less energetic 
UV photons necessary to excite the [01] and [CII], as the length 
of the material they have to pass through increases. 
Finally, there is yet a third component of material entrained in 
the outflow: that of the dust detected through the polarization 
of the visible light reported by Yoshida, Kawabata and Ohyama 
(120111) . This dust, which scatters the nuclear visible light, is 
made by grains larger and hence heavier than PAHs. In fact 
their velocity is lower (-10 km s" 1 at 1 kpc distance above the 
galaxy's disk) than that observed for the PDR component asso- 
ciated with the cold gas and presumably with PAHs. 

4.4.2. What powers the outflow in M82? 

In this Section we try to establish whether the radiation pressure 
or/and the mechanical energy due to the supernovae explosion is 
powering the PDR component participating to the motion in the 
outflow of M82 and whether the outflow is energy or momentum 
driven (i.e. whether they conserve the energy or the momentum). 

The momentum rate provided by supernovae explosion can 
be calculat ed follo wing equation 10 from Murray, Quataert and 
Thompson (120051) : 



SN 



2x 10 



33 



SFR 



1 M© yr~ 



g cm s 



-2 



(6) 



This assumes a supernovae rate vsn = 0-01 yr -1 for a 
SFR equal to 1 M© yr -1 . Forster-Schreiber et al. (2003) have 
calculated that the supernovae rate in M82 is 0.13 yr -1 and that 
the SFRs of the two starburst episodes at their peak are equal 
to 18.5 and 6.3 M© yr" 1 giving a total SFR of 24.8 M© yr" 1 . 
A supernovae rate of 0.01 yr -1 for a SFR equal to 1 M© yr" 1 
corresponds to a supernovae rate of ~ 0.25 per a SFR of 24.8 
M© yr -1 , which is a factor of 2 higher than what found in 
M82. This means that we have to divide equation 6 by 2. The 
total momentum rate available from the SN explosion is then 
~ 2.5 x 10 34 gems" 2 . 

The momentum rate given by the radiation pressure is Lsb/c, 
assuming that the optical depth at the base of the outflow 
is 1. The total bolometric starburst luminosity of M82 is 
6.6 x 10 10 L© (Forster-Schreiber et al. 2003 ), which gives a 
radiation pressure momentum rate equal to 8.5 x 10 33 g cm s~ 2 , 
a factor of 3 less than mechanical momentum rate. 
The observed momentum rate for both cones of the outflow 
(M© x Vdeproj g cm s" 2 ) ranges from 5 to 12 xlO 33 g cm s~ 2 , 
depending on the input data set (see Table [2]). 

The energy rate due to supernovae explosion under the same 
assumptions o f equa tion 6 is (equation 34 of Murray, Quataert 
and Thompson 2005]): 



E SN ~ 3 x 10 40 * 



SFR 



1 M© yr~ 



g cm s 



an 



(7) 



energy rate 



which for M82 translates in 
~ 3.7 x 10 41 ergs" 1 . 
The observed energy rate carried by the atomic PDR com- 
ponent in the outflow of M82 (\Mv 2 depro] ) ranges from 2 -5 



xlO 40 



erg s 



The momentum and energy rates in the outflow carried by 
the cold molecular component are equal to 21 x 10 33 g cm s~ 2 



and 10 41 erg s" 1 respectively, assuming a deprojected molecular 
velocity of 100 km s" 1 and a dynamical time equal to 10 7 yr 
(Walter, WeiB and Scoville 2002). We have seen in the pre- 
vious section that the cold molecular component traced by the 
CO emission is likely to be dynamically coupled to the neutral 
atomic gas PDR component. Under this assumption we can add 
the energetics of these two components to estimate a total bal- 
ance between the momentum (energy) available and those ob- 
served to move this material in the outflow. We obtained a to- 



tal observed momentum rate P a 



+ P, 



molecular 



26 - 33 x 



10 33 g cm s 2 and an observed energy deposition rate equal to 
Eatomic + E mo i ecu i ar - 12 - 15 x 10 40 erg s" 1 to be compared with 
the total momentum and mechanical energy deposition rates 
available of ~ 34 x 10 33 g cm s~ 2 and ~ 37 x 10 40 erg s" 1 re- 
spectively. 

Table |4] summarizes these values: on the left it lists the momen- 
tum (energy) rate available in the system and on the right the ob- 
served momentum (energy) rates in M82. Both the total momen- 
tum and energy deposition rates available in the system are com- 
parable to the observed values. The momentum driven outflow 
is not compatible with the observations if the source of power 
is the radiation pressure alone. The energy driven case seems to 
leave more room for agreement between the observed and the 
available energy of the PDR outflowing material. 
However, if we include in the energy and momentum budget 
calculation also the numbers necessary to sustain the outflow- 
ing material observed in Ha by Shopbell and Bland-Hawthorn 
(1998) the energy driven case becomes only marginally com- 
patible with the observations. In fact, assuming a mass of the 
outflow observed at these wavelengths equal to 5.8 x 10 6 M©, 
a kinetic energy equal to 2.0 x 10 55 erg, a mean deprojected ve- 
locity of ~ 600 km s" 1 ( Shopbell and Bland-Hawthorn [T9981) 
and a length of the outflow above the galactic plane equal to 1 
kpc, one derives a momentum rate ~ 2.2 x 10 33 g cm s~ 2 and an 
energy rate ~ 40 x 10 40 erg s" 1 (also listed in Table [4}. Taking 
into account the overall uncertainties (~ a factor of 2) of the 
numbers derived from the above calculation and considering the 
cold molecular gas, the neutral atomic gas and hot ionized gas 
phases energetic together, we conclude that, assuming a simple 
jet-like outflow geometry, both the energy and momentum driven 
process are compatible with the observed values, although the 
momentum driven case seems to be slightly favored. 



4.5. Line diagnostic diagrams 



Figure 21 shows the ratio of each of the observed lines to the FIR 
continuum as well as the [OI] 63 yum/[CII], and [OIII]/[OI] 63 
yum ratios as a function of the F^m]/ ^[i(%m] ratio, equivalent to 
the IRAS colors, derived from our PACS maps. The latter traces 
the dust temperature and therefore, to first order also the inter- 
stellar radiation field. The M82 values are shown in gray scale. 
Overplotted are the data collected from the literature or mea- 
sure d in th e PACS GT Key Program SHINING (Gracia Carpio 
et al. 1201 lb . Different symbols refer to different galaxy types (see 
figure caption for details). The emission line data points of M82 
are those derived from the final maps reduced to the same spa- 
tial resolutions. The maps at 60 and 100 yum have been derived 
as explained in Section 4.1.1. 

There is excellent agreement between the M82 data, which probe 
regions of linear extent of ~ 300 pc, and the global emission 
properties of galaxies. This means that the physical mechanisms 
governing the observed relationships act up to a scale as small as 
300 pc, and they are the same independent on whether they act 
in the starburst region, the diffuse disk or the outflow. The only 
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Table 4. Available and observed molecular and atomic gas energetics in the outflow of M82 



Momentum Rates 


Available 




Used 




10 33 g cm s -2 




10 33 g cm s -2 


p 

*- S N 


25 


p 

' molecular 


21 




8.5 


p 

1 atomic 


5-12 






P Ha 


2 


total 


33.5 




28-35 


Energy Rates 


Available 




Used 




10 40 g s" 1 




10 40 g s" 1 


Esn 


37 


Emoeluclar 


10 






F 

atomic 


2-5 






Ehcx 


40 


total 


37 




52-55 



difference is in the panel showing the [01] 63 pm /FIR ratio that 
in the case of M82 is remarkably constant, while the galaxy sam- 
ple is exhibiting a large scatter. It is interesting to note that the 
scale we probe (~ 300 pc) is comparable to the upper limit of the 
scale where the Kennicutt-Smidth law begins to break (Schruba 

etaLnnni. 

Each left and middle panel of Figure [22] shows the 
FIRu ne IFIR continuum ratios versus the F [60 ^m]/^[i(%m] colors and 
the FIR luminosities for M82 only, color coded depending on 
whether the regions belong to the outflow, the starburst or the 
disk of M82 as shown in the upper panel of Figure 18 The top 
panels show the [CII]/FIR ratio decreases by one order of mag- 
nitude with the increase of the F^O/im]/ ^[loo/zm] IRAS color and 
the FIR luminosity, nicely following the decreasing trend out- 
lined from the whole SHINING galaxy sample (Gracia Carpio 
et al. [20TTb . As expected, the lowest [CII]/FIR values in M82 
arise from the starburst region (green points). This decrease was 
first discovered with the ISO LWS and it has been referred to as 
the [CII]-defic it (Fischer et alg999l[20 QTl Ma lohtra et al. l200Tl 
Luhman et al. [T9981 [20031 Abel et al. [20091) . i.e. the warmer, 
more active and more infrared luminous the systems are, the less 
[CII] with respect to their FIR emission they emit. 

The near universality of the [CII] deficit for luminous 
systems has been recently challenged by observations of high 
redshift galaxies, which show a "normal" [CII]/FIR ratio 
despite of their high FIR luminosities (Maiolino et al. 2009, 
Haile y-Dunsheath et al l20T0l Iviso n et al . l20T0l Stacey et al. 
l20T0l ,Wagg et al. l20T0l Sturm et al. l20T0b . On the other hand, 
if the [CII]/FIR ratio is plotted against the ratio between the 
FIR luminosity and the H2 molecular mass, rather than FIR 
luminosities, the scatter in the relation is significantly reduced 
as is the disparity between systems at low and high redshift. 
This means that the [CII] deficit becomes universal at all 
redshifts above a threshold of the L FIR /M Hl ratio, a commonly 
used indicator of star formation efficiency or evolutionary stage 
(Gonzalez-Alfonso et al, [20081 Gracia Carpio et al, 1201 11 ). 
Furthermore, Gracia Carpio et al. find that this deficit extends 
to other PDR and pure HII region FIR lines (see their figure 
2). In galaxies where L FIR is predominantly powered by star 
formation, the L FIR /M Hl ratio is, to first approximation, the 
ratio between the energy released by the star formation and 
the gas reservoir from where the stars form and therefore it is 
a parameter directly related to the star formation efficiency. 
This is not the case for the L FIR . Gracia Carpio et al. (120111) 
show that above a certain L FIR /M H2 threshold one also finds an 
increase of the ionization parameter U, i.e. the ratio between 
the ionizing photons to the gas density at the inner surface 
of the ionized cloud. The higher this parameter is the larger 



the HII region and the FIR emitted by the dust heated in 
the HII region. This increasing FIR does not mean a corre- 
sponding increase in the FIR emission lines arising from the 
PDR, thus the global FIR(PDR) line / FIR(HII + PDR) continuum 
ratio decreases. The values of the L FIR /M Hl threshold (> 80 
L Q M~ l ) is very close to the threshold above which the star 
formation becomes much more efficient than what is predicted 
by the Kennicutt- Schmidt relation for normal star forming 
galaxies (Genzel et al. 12010b . This similarity suggests that 
the galaxies and the regions within galaxies which exhibit a 
FIRu ne deficit, are in this highly efficient mode of star formation. 

In order to determine in which mode of star formation the 
various regions of M82 are, we want to plot the line/FIR ratios 
of the points belonging to M82 as a function of L FIR /M H2 . Since 
we do not have a CO map of M82 available from which we can 
derive the H2 mass, we use instead the mass of the PDR gas 
to calculate the L FIR /M[cii](H) ratio. The cold molecular gas is 
more closely related to the site of star formation than the atomic 
neutral gas, but the atomic gas we are considering here is that as- 
sociated with the PDR and therefore it is also related to the star 
formation, although not as directly as the molecular gas. This is 
why the ratio between the FIR luminosity and the mass of neu- 
tral hydrogen associated with the PDRs can be considered to first 
order, as star formation efficiency. 

We have produced plots analogous to the FIRu ne / FIR continuum 
deficit diagrams shown in Gracia Carpio et al. (2011) for M82 
(color coded: disk in black, starburst in green, north outflow in 
cyan and south outflow in red, see also the schematic figure on 
top right of Figure 18). These diagrams are shown in the right 
panels of Figures |22|for each line. We have also plotted the hori- 
zontal lines that Gracia Carpio et al. ( 201 1 ) consider as the limit 
between galaxy with normal and deficient F IRu ne I F IR continuum 
ratios. 

The first thing to notice is that even the starburst region 
of M82 is not in this high efficiency mode of star formation, 
since all points are above the horizontal lines. Nevertheless we 
observe different behaviors among the different lines. 
The [CII]/FIR ratio gently decreases in all cases, even when 
plotted as function of the L FIR /M[ C ii](H), whereas the similar 
plot in Gracia Carpio et al., (2011) who used the mass of H 2 , 
shows an almost constant behavior up to (~ 80 L© M~ l ), and 
then a decrease. However, we notice that the data shown in 
Gracia Carpio (I20TT1) did not include the whole SHINING 
sampl^] and the data available at that time contained more 



3 The SHINING sample, P.L E. Sturm, includes 50 galaxies ranging 
from HII, starbursts, Seyfertl and Seyfert2, and LINERs and including 
LIRGs and ULIRGs, Gracia Carpio et al. (in preparation) 
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ULIRGs/AGNs than starbursts and star forming galaxies. Now 
that the sample is complete and more starbursts/star forming 
galaxies populate the low L FIR /M Hl part of the plot, a decrease 
at low L fir /Mh 2 is also clearly visible for the SHINING sample 
(Gracia Carpio et al. in preparation) although with a slope much 
less steep than that followed by galaxies at L FIR /M H2 > 80. This 
continuous decrease is predicted by models which assume an 
increasing production of FIR emission arising by the dust heated 
in the HII regions, as the ionization parameter U increases (see 
Figure 12 of Abel et al. [20091) . 

The [01] 63 jim /FIR and the [01] 145 /im /FIR ratios in 
M82 are remarkably constant and, in contrast to the [CII]/FIR 
behavior, do not show any hint of decline in M82. The constancy 
is compatible with what is observed in the whole SHINING 
sample, which shows a decline of the [OI]/FIR ratio only at 
60/100 ratios higher than that probed in M82. 

Both trends can be understood by modeling these ratios 
with CLOUDY as shown in Figure 3 of Gracia Carpio et al. 
(2011). These figures show the FIRu ne /FIR ratios on diagrams 
of the ionization parameter U versus the gas density n H . Each 
diagram is divided into gray and red zones, corresponding to the 
typical values of starburst/star forming galaxies and ULIRGs 
respectively. The loci of the line over FIR value corresponding 
to the transition between the zones correspond to the horizon- 
tal lines shown in Figure 22 First we address the diagrams 
involving only the PDR lines ( [CII], and the two [OI] lines): 
when moving diagonally on these diagram, i.e. increasing both 
parameters, the [CII]/FIR ratio decreases but the [OI]/FIR ratios 
do not change significantly because in the first case one moves 
perpendicular to the loci of constant [CII]/FIR and in the second 
case one moves almost parallel to the loci of constant [OI]/FIR, 
at least in the range of line/FIR ratios spanned by M82 and 
starbursts in general (gray zones in Figure 3 of Gracia Carpio 
et al. 1201 II) . The behavior of the lines arising from ionized gas 
is reversed with respect to what we have just examined in the 
case of PDR lines. In particular we concentrate here on the 
[OIII]/FIR line loci, on the U-nu diagram. We can see that 
in this case, increasing both parameters means increasing the 
HII-line/FIR ratio, as it is observed in M82 shown in the bottom 



panel of figure 22 However, it is worth mentioning that the pure 
HII region lines are very sensitive to the input SED used for 
the modeling, and in particular to the hardness of the exciting 
photons, much more than the lines originating from PDRs or the 
dust temperature and total flux which is sensitive to the energy 
of the heating photons rather than to the hardness. 
From the U-nn diagram we can derive the FIRu ne / FIR continuum 
ratios expected for the values calculated by Forster Schreiber et 
al.. l200T1 for M82 (U = -2.3 dex and n e ~ n H = 300 cm" 3 ). 
These values agree well to the averaged values we observe, as 



shown by the green crosses in all panels of 22 



Finally, a closer look to the trend followed by the 



[Oil I] I FIR ratio in Figure 22 seems to suggest that there ex- 
ist two different behaviors in M82. In the disk (black points) 
and in the starburst (green points) this ratio increases with IRAS 
color, FIR luminosity and L FIR /M Hl as expected, (see inner pan- 
els in the bottom row in Figure |22| ). On the other hand, in the 
outflow this ratio remains constant and is higher than the corre- 
sponding point of the disk at low IRAS color, FIR emission and 
L FIR /M H2 ratios. This suggests that in the outflow the [OIII] is 
enhanced with respect to the FIR emission. We have seen that 
the [OIII] emission in the outflow is dynamically coupled with 
the gas emitting in the [CII] and [OI] lines, that we have asso- 



ciated to PDR clouds. We have also seen that the PDRs in the 
outflow survive in form of small cloudlets. We can imagine that 
the dust associated with these small clouds is shielded from the 
UV photons resulting in colder FIR colors and lower FIR emis- 
sion. We notice that shocks do not play a significant role in the 
heating of the PDR material anywhere in M82 including the out- 
flow, since the [OI] 63 /mi /[CII] ratio is always well below the 
values (~ 10) predicted in the presence of shocks (Hollenbach 
and McKeeQSHI]). 

5. Summary and Conclusions 

We have mapped a 2.5' x2. 5 ' region of M82 in four FIR atomic 
fine- structure lines, [CII] at 158 yum, [OI] at 63 and 145 //m 
and [OIII] at 88 /mi, with the PACS spectrometer on board the 
Herschel Space Observatory, reaching unprecedented spatial res- 
olution (~ 200 - 300 pc) at these wavelengths. The mapped area 
covers approximately 1 kpc along the minor axis of the galaxy. 
This enables detection of the M82 outflow structures in these 
lines for the first time. Since the [CII] and [OI] lines are the main 
coolants of the cold neutral medium, these observations trace 
and probe the atomic gas in the outflow. This paper presents the 
first analysis of these observations with emphasis on the nature 
and structure of the atomic gas in the outflow and on its relation 
with other ISM phases studied in previous works. We summarize 
the main results as follows. 

- We derived intensity, radial velocity and dispersion maps of 
the entire mapped region in all lines. We also used the offline 
continuum images and the parallel array scans to produce a 
dust temperature map. The line dispersion maps most clearly 
show two dynamically separated structures, disk and out- 
flow, since the dispersion in the outflow is ~ 50 - 100 km s" 1 
higher than that in the disk. These distinct components are 
also visible in the radial velocity maps and indicate that 
the northern outflow is receding from us while the south- 
ern outflow is moving toward us, in agreement with previ- 
ous results. In the outflow, the mean deprojected velocity, 
calculated assuming a simply inclined jet-like geometry, is 
~ 75 km s" 1 for both the cold neutral and the ionized gas 
traced by the [OIII] 88 jim line, close to the average veloc- 
ity of the molecular gas (~ 100 km s" 1 ) reported by Walter, 
WeiB and Scoville (2002) and much smaller than the veloc- 
ity of the outflowing material observed in Ha (~ 600 km s" 1 , 
Shopbell and Bland-Hawthorn [1998]) 

- Analysis of the line ratio maps and of the line-to- continuum 
maps confirms earlier findings that the base of the northern 
outflow is more obscured than the southern one due to the 
inclination of the galaxy disk. These maps also reveal that 
the opening angle of the cold gas component of the outflow 
is much larger than the ionized gas traced by the [OIII] 88 
yum which appears much more collimated. 

- We have carried out PDR modeling of the neutral gas emis- 
sion lines to derive the key physical parameters of the atomic 
gas within and surrounding the starburst, stellar disk and the 
outflow components of M82. By applying the Kaufman et 
al. (1999 ) PDR models to each pixel of our line intensity 
images we were able to produce, for the first time in an ex- 
ternal galaxy, maps of the FUV interstellar radiation field 
Go, the gas density n H , the gas temperature T gas , the atomic 
gas cloud beam filling factor, the column density of atomic 
hydrogen, and the [CII] opacity. Since the correction for the 
contribution of [CII] emission from ionized gas is very un- 
certain we run the model with and without it, but the re- 
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Fig. 21. [CII]/FIR, [01] 145yum/FIR, [OIII]/FIR, [01] 63 yum/FIR, [01] 63 yum/[CII] and [OIII] 88/[OI] 63 //inversus the continuum 
60/100 jim ratio. M82 values are shown as gray scale. Symbols are the following (Gracia Carpio et al. 1201 H and in preparation): Big 
symbols: SHINING/PACS observations; Filled symbols: the lines are detected; Open symbols: one or both lines are not detected; 
Blue points: HII galaxies; Green points: LINER galaxies; Black symbols: AGNs; Circles: Seyfert 2 galaxies; Triangles: Seyfert 1 
galaxies; Diamonds: QSOs; Magenta points: high-z galaxies; Grey points: unclassified galaxies; Small magenta points: flux at 25 
or at 100 yum is an upper limit. 



suiting maps do not differ significantly. Solution maps using 
each [OI] line together with the [CII] and continuum TIR 
emission were similar to within the errors for all parame- 
ters except the gas density, probably because of alignment 
problems between the [01] 63 /mi and [CII] map. For this 
reason and also because the [01] 63 jim optical depth can 
vary significantly in the galaxy we based our analysis only 
on the results obtained using as input to the model the [01] 
line at 145 jim together with [CII] and TIR. The area beam 
filling factor decreases by up to an order of magnitude in the 
diffuse disk and in the outflow. Go and hh both vary from 



~ 50 to few 10 3 and the atomic gas temperatures are in the 
range of 300-400 ° K. The [CII] optical depth is everywhere 
much less than unity, enhanced along the disk but decreasing 
toward the starburst region, perhaps because some material 
has been cleared out by the winds that power the outflow. 
- The relation between Go and n H in M82 is the same as has 
been observed in previous works for the global emission of 
star forming galaxies extended by almost two orders of mag- 
nitude towards smaller values. In this plot, the same relation 
is followed regardless of whether the points belong to the 
starburst, the outflow or the diffuse disk, confirming that the 
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Fig. 22. The Line/FIR ratio versus: left: 60/100 colors; middle: L FIR ; right : L FIR /M(HI) PDRs . Data are for M82 only. Black points 
are from the diffuse disk; red points from the southern outflow; cyan points from the northern outflow and green points from the 
starburst. The small frames in the last row show the same relation reported in the big frames but only for the points belonging to 
the disk (black) and the starburst (green). The horizontal lines correspond to the threshold identified by Gratia Carpio et al. (2011) 
between galaxy with normal (above) and deficient (below) FIRi ine /FIR continuum ratios, or between normal and high efficency of SF. 
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clouds in the outflow are still organized in PDRs very similar 
to the clouds in the rest of the galaxy. 
We find that the mass and the kinetic energy of the outflow- 
ing cold neutral gas are > 2 - 8 x 10 7 M© and > 1 - 5 x 10 54 
erg slightly smaller than the corresponding values found for 
the outflowing molecular gas (~ 3.3 x 10 8 M©, ~ 3 x 10 55 
erg). The mass outflow rate of the molecular and atomic gas 
is ~ 43-58 M yr _1 comparable to the SFR of M82 equal to 
25 M yr _1 , resulting in a mass loading factor M/SFR ~ 2. 
This result, together with the fact that these two components 
show similar velocities, suggests that they are dynamically 
coupled with similar origin. Since the cold atomic gas is 
consistent with emission from classical PDRs (i.e. interface 
regions between the ionized and molecular gas) with small 
filling factors, the coupling of these two cold components 
of the ISM is naturally explained if both the molecular and 
atomic media belong to the same disk clouds entrained in 
the outflow by the wind where they partially evaporate due 
to heating driven by thermal conduction from the hot gas, 
surviving as clouds smaller than their original size. This sce- 
nario is supported by the fact that the decline of the derived 
FUV ISRF in the outflow is consistent with a pure geomet- 
rical dilution, indicating that FUV photons can travel in the 
outflow with a low probability of being absorbed by dust, 
which therefore must be concentrated in small clouds. This 
also implies that in the outflow, there is no need for a signif- 
icant in situ star formation in addition to the starbust light. 
We have found that in the outflow, the ionized gas traced by 
the FIR [OIII] 88 /mi line flows at a much smaller velocity 
(~ 75 km s" 1 ) than the ionized gas traced by the Ha emission 
(~ 600 km s"\ Shopbell and B land-Hawthorn [19981) but at a 
velocity similar to that of the neutral atomic gas traced by the 
[CII] and [01] lines. We have also not found line splitting in 
any of the observed FIR lines, including the [OIII] line at 88 
/mi, despite the fact that the line splitting detected in Ha is 
large enough to be detected at the PACS spectral resolution. 
These results suggest that [OIII] 88 /mi and Ha line emis- 
sion trace two different ionized gas components, both in ve- 
locity and in physical space. We propose a scenario in which 
the bulk of the [OIII] line at 88 /mi traces gas surrounding 
the PDR cloudlets observed in the [CII] and [01] lines, pho- 
toionized by the starburst radiation. This would justify the 
fact that all these components flow at the same velocity. In 
this scenario, the bulk of the Ha emission arises from the 
shocked gas at the interface between the hot plasma wind 
emitting in X-ray and the halo. This would explain why the 
Ha emission is mainly located in the walls of the outflow as 
demonstrated by the observed line splitting. 
The momentum and the energy available to drive the out- 
flowing gas, molecular atomic and ionized, are both compat- 
ible with that observed in the flows, although the momentum 
driven case seems to be slightly favored by the data. 
We compare the trends found in the relative FIR fine- 
structure line fluxes as a function of the IRAS 60/100 col- 
ors of normal and infrared-bright galaxies to those found in 
the spatially resolved, ~ 300 pc regions of M82. We find 
that for all but the warmest far-infrared colors of ULIRGs 
and the most extreme LIRGs, the ratios in M82 overlap with 
the trends found recently by Gracia Carpio et al. (120111) for 
the whole galaxies. These authors explore parameter space 
in density and ionization parameter. Based on photoioniza- 
tion models, they argue that the atomic fine- structure line 
deficits found in the neutral gas of galaxies with the warmest 
FIR colors and the highest values of L FIR /M Hl are caused 



by the presence of high ionization parameters in which the 
dust rather than the gas absorbs most of the ionizing pho- 
tons. They suggest that these most extreme galaxies are in 
the high efficiency mode of star formation rather than in a 
normal/Milky Way-like mode of star formation. Importantly 
we find that down to 300 pc size regions, even in regions of 
the highest star formation efficiency in M82, only the normal 
mode of star formation is present, i.e. no hint of bimodality 
is found. 
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Appendix A: 

We have seen in Section 4.2 that we find a significant difference 
in the morphology of the gas density map obtained using as in- 
put to the PDR modeling the datasets with the two [OI] lines 
(Figure [15]), and this difference is independent on the correction 
we apply to the observed [CII] emission due to the contribution 
from ionized gas. If this difference were caused by the fact that 
[01] at 63 jim is thick in these regions, we would expect to ob- 
serve less [01] 63 than what is predicted by the solution obtained 
using the [01] line at 145 //m, which is optically thin. We have 
checked this and we find the opposite answer, that is in order to 
make the two density map solutions similar, we should reduce 
the input flux at 63 /mi. In other words, it looks like, in these re- 
gions we have too much [01] 63 /mi flux. Although this could be 
related to the uncertainties of the geometrical factor equal to two 
we have applied to the observed [01] 63 line flux (see section 
4.1.2), we wanted to check whether this difference could arise 
from our processing steps. 

Each set of input values to the PDR modeling were obtained 
from maps (either [CII], FIR and [01] at 63 /mi or [CII], FIR and 
[01] at 145 /mi), rebinned to the largest pixel size (6", that in the 
red channel), smoothed to the worse resolution (that of the [CII] 
map), reduced to the same size and aligned to the [CII] peak. In 
this chain of processing, the map which suffers the heaviest ma- 
nipulation is that at 63 /mi. This is because the original pixel size 
is smaller, and its PSF differs from that at 158 /mi much more 
than that at 145 /mi. Moreover, the observations at 158 and 145 
/mi were carried out in the same AOR, therefore the pointing ac- 
curacy is likely to affect these two observations in a similar way. 
This may give rise to alignment problems between the final map 
at 63 /im and those at higher wavelengths. 
For these reasons we decided to run the models again, this time 
using the [01] 63 /mi map aligned on the [CII] map on the steep 
gradient of the flux distr ibution rather than to the peak (shift of 1 
pixel = 6"). Figure |AJ"| clearly shows that now both gas density 
maps show an enhancement toward the north part of the star- 
burst, although the values are lower for the solutions obtained 
with the [01] 145 //m line than those obtained with the [01] 63 
/mi line. The maps in the other parameters, Go, T Ras and 0, show 
the same morphology as those shown in figures 14 16 and 17 
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Fig. A.l. The PDR solution maps obtained for the gas density uh using the [01] map at 63 \i shifted to match the steep gradient of 
the starburst region to that in the [CII] map, instead of the [CII] peak emission. 
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